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> ■ ABSTRACT 
(N 

^ . We have discovered 16 Type la supernovae (SNe la) with the Hubble Space 

O . Telescope (HST) and have used them to provide the first conclusive evidence 

Q . for cosmic deceleration that preceded the current epoch of cosmic acceleration. 

^ . These objects, discovered during the course of the GOODS ACS Treasury 

o 

^ . 

c/3 . distances to these objects, and to 170 previously reported SNe la, have been 

^ ' determined using empirical relations between light-curve shape and luminosity. 



program, include 6 of the 7 highest-redshift SNe la known, all at z > 1.25, 
and populate the Hubble diagram in unexplored territory. The luminosity 



A purely kinematic interpretation of the SN la sample provides evidence at 
the > 99% confidence level for a transition from deceleration to acceleration 
or similarly, strong evidence for a cosmic jerk. Using a simple model of the 
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expansion history, the transition between the two epochs is constrained to be at 
z = 0.46 ± 0.13. The data are consistent with the cosmic concordance model 
of Qm ~ 0.3, f^A ~ 0.7 {Xdof ~ 1-06), and arc inconsistent with a simple model 
of evolution or dust as an alternative to dark energy. For a fiat Universe with 
a cosmological constant, we measure Qm = 0.29±qq3 (equivalently, Q\ = 0.71). 
When combined with external fiat- Universe constraints including the cosmic 
microwave background and large-scale structure, we find w — — 1.02±q ;Jg (and 
w < —0.76 at the 95% confidence level) for an assumed static equation of state 
of dark energy, P — wpc^. Joint constraints on both the recent equation of 
state of dark energy, wq, and its time evolution, dw/dz, are a factor of ~ 8 
more precise than its first estimate and twice as precise as those without the 
SNc la discovered with HST. Our constraints are consistent with the static 
nature of and value of w expected for a cosmological constant (i.e., wq = —1.0, 
dw/dz = 0), and are inconsistent with very rapid evolution of dark energy. We 
address consequences of evolving dark energy for the fate of the Universe. 

subject headings: galaxies: distances and redshifts — cosmology: observations — 
cosmology: distance scale — supernovae: general 

1. Introduction 

Observations of type la supernovae (SNe la) at redshift z < 1 provide startling and 
puzzling evidence that the expansion of the Universe at the present time appears to be 
accelerating, behavior attributed to "dark energy" with negative pressure (Riess et al. 1998; 
Perlmutter et al. 1999; for reviews, see Riess 2000; Filippcnko 2001, 2004; Leibundgut 
2001). Direct evidence comes from the apparent faintness of SNe la at 2; O.5.. Recently 
expanded samples of SNe la have reinforced the statistical significance of this result (Knop 
et al. 2003) while others have also extended the SN la sample to z ~ 1 (Tonry et al. 
2003; Barris et al. 2004). Observations of large-scale structure (LSS), when combined with 
measurements of the characteristic angular size of fluctuations in the cosmic microwave 
background (CMB), provide independent (though indirect) evidence for a dark-energy 
component (e.g., Spergel et al. 2003). An independent, albeit more tentative investigation 
via the integrated Sachs- Wolfe (ISW) effect also provides evidence for dark energy (Scranton 
et al. 2003). The magnitude of the observed acceleration was not anticipated by theory 
and continues to defy a post facto explanation. Candidates for the dark energy include 
Einstein's cosmological constant A (with a phenomenally small value), evolving scalar 
fields (modern cousins of the inflation fleld; Caldwell, Dave, & Steinhardt 1998; Peebles 
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& Ratra 2002), and a weakening of gravity in our 3 + 1 dimensions by leaking into the 
higher dimensions required in string theories (Deffayet, Dvah, & Gabadadze 2002). These 
explanations bear so greatly on fundamental physics that observers have been stimulated 
to make extraordinary efforts to confirm the initial results on dark energy, test possible 
sources of error, and extend our empirical knowledge of this newly discovered component of 
the Universe. 

Astrophysical effects could imitate the direct evidence from SNe la for an accelerating 
Universe. A pervasive screen of grey dust could dim SNe la with little telltale reddening 
(Aguirre 1999a,b). Luminosity evolution could corrupt the measurements if SNe la at 
z ~ 0.5 are intrinsically fainter than their low-redshift counterparts. To date, no evidence 
for an astrophysical origin of the apparent faintness of SNe la has been found (Riess 2000; 
Coil et al. 2001; Leibundgut 2001; Sullivan et al. 2003). However, given the significance 
of the putative dark energy and the unique ability of SNe la to illuminate it, we need a 
more definitive test of the hypothesis that supernovae at z ~ 0.5 are intrinsically dimmer, 
or dimmed by absorption. 

If cosmic acceleration is the reason why SNe la are dimmer at z ~ 0.5, then we expect 
cosmic deceleration at 2; > 1 to reverse the sign of the observed effect. The combination of 
recent acceleration and past deceleration is a clear signature of a mixed dark-matter and 
dark-energy Universe and one which is readily distinguishable from simple astrophysical 
dimming (Filippenko & Riess 2001). 

Furthermore, assuming SNe la at 2; > 1 continue to trace the cosmological world 
model, measurements of SNe la in the next redshift octave provide the unique ability to 
discriminate between a static and evolving dark-energy equation of state. This would 
provide a vital clue to distinguish a cosmological constant from other forms of dark energy 
that change with time. 

Ground-based efforts to look for past deceleration with SNe la have offered hints of 
the effect, but ultimately they have suffered from insufficient signal-to- noise ratios (Tonry 
et al. 2003; Barris et al. 2004). Discovering, confirming, and then monitoring transients 
at l ^ 25 mag on the bright sky is challenging even with the largest telescopes and the 
best conditions. A single SN la at z 1.7, SN 1997ff, discovered with WFPG2 on the 
Hubble Space Telescope (HST) (GiUiland, Nugent, & Phillips 1999), provided a hint of past 
deceleration; however, inferences drawn from a single SN la, while plausible, are not robust 
(Riess et al. 2001; Benitez et al. 2002; Mortsell, Gunnarsson, & Goobar 2001). 

To study the early expansion history of the Universe, we initiated the first systematic, 
space-based search and follow-up effort to collect SNe la at 2; > 1, carried out in conjunction 
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with the Great Observatories Origins Deep Survey (GOODS) Treasury program (Giavalisco 
et al. 2003) conducted with the Advanced Camera for Surveys (ACS) aboard HST. (The 
abihty to detect SNe at z > 1 with the Space Telescope was an apphcation first envisioned 
during its planning; Tammann 1977, Colgate 1979). A separate "piggyback" program was 
utilized to obtain target of opportunity (ToO) follow-up HST observations of the SNe la 
with ACS and NICMOS (the Near- Infrared Camera and Multi-Object Spectrograph). 
Elsewhere we present a color-based method for discrimination of SNe la at z > 1 from other 
transients (Riess et al. 2003) and the full harvest of the SN survey (Strolger et al. 2004). 
We present the follow-up spectroscopy and photometry of 16 SNe la in §2, hght-curve 
analysis in §3, cosmological tests and constraints in §4, and a discussion and summary in 
§5 and §6, respectively. 

2. Target of Opportunity Follow-up: Light Curves and Spectra 

2.1. Follow-up 

The methods and criteria we used to search the GOODS ACS Treasury data for SNe 
are described by Strolger et al. (2004) and are based on image subtraction (Perlmutter 
et al. 1997; Schmidt et al. 1998). Strolger et al. (2004) provide the parameters of the 
search including search depth, efficiency, timing, and false-positive discrimination, as well 
as a hst of all detected SNe. Briefly, our search was conducted in the F850LP (Z-band) 
to an effective hmit of ~26.0 (Vega) magnitude covering 0.1 square degree in 5 epochs (at 
intervals of ~45 days). Our limiting magnitude was 1 to 2 mag fainter than the expected 
peak of a SN la over the target range of 1 < 2; < 1.6, therefore SNe la we collected (whose 
intrinsic dispersion is expected to be < 0.2 mag)would not preferentially be selected from 
the bright tail of their intrinsic distribution. In Table 1 we provide discovery data for the 
SNe la reported here. 

Our ToO candidates were generally too faint to anticipate useful spectral discrimination 
from the ground; it was therefore necessary to initially identify SNe la photometrically. To 
discriminate SNe la at 2; > 1 from SNe II and from SNe I at lower redshifts, we used a 
combination of photometric redshifts of the host galaxies (with 9 passbands) and rest-frame 
ultraviolet (UV) colors; see Riess et al. (2003) for details. A comparison of the photometric 
redshifts available at the start of the survey with spectroscopic redshifts obtained thereafter 
yields an RMS of 0.05 for the quantity {zphot — -2spec)/(l + ^^^pec) for 25 hosts (with a 4 cr 
outlier; see Strolger et al. 2004). 

We selected 9 individual candidates for subsequent ToO observations (including one 
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pair of targets observed within the same ACS field). In addition to these primary targets, 
judicious positioning of the follow-up fields provided serendipitous monitoring of 4 additional 
high-redshift SN la candidates. For 4 more SNc la, the periodic imaging of the GOODS 
survey (sometimes augmented with a few ground-based observations) provided sufficient 
characterization of their light curves. Subsections of the discovery and pre-discovery images, 
as well as their difference (centered on each SN la), arc shown in Figure 1. 

The great benefit of HST observations-dramatically reduced sky noisc-is only fully 
realized if operational constraints can be overcome. To get the full advantage of HST, the 
supernova search needs to be completed, the best targets selected, spectra taken, and the 
photometric follow-up set in motion in 2 weeks. Otherwise, the peak of the light curve will 
be missed, the dechne rate of the object will not be well measured, and, as the objects 
fade, the spectra will become too difficult to obtain. Our challenge was to do all this 
without using highly disruptive and inefficient "24 hour" ToOs (for which the HST weekly 
schedule is immediately interrupted, wasting precious observing time). To achieve our 
goals, each GOODS epoch was scheduled during a 3-4 day interval immediately preceding 
the weekly deadline for a ToO activation. This allowed the ToO observing to be scheduled 
and uploaded to the spacecraft for the following week's observations without delay. By 
adopting this prescription we were able to achieve a time interval of typically 9 to 11 days 
between the discovery and first follow-up observation, a span of less than 5 days in the rest 
frame of a SN at 2; ;^ 1. 

After discovery and photometric screening, we developed a follow-up plan which was 
appropriate for our best estimate of the SN redshift. AGS and the F850LP filter were 
typically used for 6-8 epochs to obtain a rest- frame S-band or [/-band light curve for 
SNe la at 2; > 1 extending to ~ 20 rest-frame days past maximum brightness. This strategy 
utilized results from Jha, Riess, & Kirshner (2004a), who demonstrated the utihty of a large 
sample of [/-band light curves for constraining the light-curve shape parameter and phasing 
of the optical light curves. To provide rest-frame optical zero-points, where SNe la are best 
calibrated, we used NIGMOS Gamera 2. SNe expected to be at 1.0 < 2; < 1.3 were imaged 
with NIGMOS Gamera 2 and FllOW. SNe expected to be at 2; > 1.3 were imaged with 
both FllOW^ and F160W. The one exception was SN 2002ki, for which a ground-based 
spectrum yielded z = 1.14 from narrow, host-galaxy [O 11] emission, and F160W was used 
to sample the second maximum in the rest-frame infrared (IR) to use as a secondary means 
of classification. 

AGS coupled with a grism filter (G800L) provides slitless spectroscopy over the entire 
wide-field camera (WFC) field-of-view. For obtaining spectra of high-redshift SNe la, 
this mode of observing has noteworthy advantages and disadvantages as compared to 
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spatial-blocking spectroscopy with STIS on HST or from the ground with a large-aperture 
telescope. The primary advantage of ACS grism spectroscopy is its efficiency. ACS is 
the most efficient camera to fly on HST due to its rationed, silver-coated reflections. 
Importantly, slitless spectroscopy with the grism retains the efficiency of HST resolution 
along one spatial direction, dispersing the light of a point-spread function (PSF) over few 
sky pixels. Another advantage attuned to identifying high-redshift SNe la is the relatively 
low sky brightness between 0.8 /xm and 1 fim from the vantage point of HST. 

Using the HST grism has disadvantages as well. HST is a factor of 3-4 smaller in size 
than the largest ground-based optical telescopes. Moreover, the grism disperses a large 
spectral range of the sky onto the position of the SN. For G800L, wavelengths shorter than 
~5500 A are blocked so the total sky counts are only ~50% greater (hence ~25% more 
noise) than for the /-band (F814VF). A more troublesome feature of grism observations is 
the superposition of multiple sources from multiple diffraction orders on the same sky pixels. 
Careful consideration must be given to possible contamination of a SN spectrum by nearby 
sources, especially the host galaxy (although this effect can be mitigated by judicious choice 
of the telescope roll angle). The ACS grism also is limited in resolution to R = A/AA = 200. 
However, this resolution is well-matched to measure SNe la whose blended absorption 
features are broadened by ejecta dispersions of lO'' km s~^. Serendipitous spectra of SNe la 
obtained during ACS commissioning and reported by Blakeslee et al. (2003) provided the 
first in-orbit examples at high redshift (specifically, SN 2002dc at z = 0.47 and SN 2002dd 
at ^ = 0.95). 

We employed the ACS grism for our ToOs when we expected to obtain a sufficient 
signal-to-noise ratio (S/N) for classification and redshift determination in fewer than 8 orbits 
of integration time and without detrimental contamination. These included SNe 2002fw, 
2003az, 2003dy, 2003es, and 2003eq as well as the neighboring SN 2003eb. Spectra of other 
SNe or host galaxies were obtained with ground-based telescopes (Table 3), including the 
VLT, Magellan, Keck-II with NIRSPEC (McLean et al. 1998), and especially Keck-I with 
LRIS (Oke et al. 1995). 



2.2. Photometry 

After the search phase, all images were reprocessed using up-to-date reference files 
and the CALACS pipeline in the STSDAS package in IRAF^. This procedure includes 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
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Table 1. Discovery Data 



SN 


Nickname 


UT Date 


SN a(J2000) 


SN ^(J2000) 





2002fw 


Aphrodite 


Sep. 19.86 


03:32:37.52 


-27:46:46.6 


2002£x 


Athena 


Sep. 20.84 


03:32:06.80 


-27:44:34.4 


2002hp 


Thoth 


Nov. 1.51 


03:32:24.79 


-27:46:17.8 


2002hr 


Isis 


Nov. 1.64 


03:32:22.57 


-27:41:52.2 


2002kc 


Bilbo 


Dec. 21.50 


03:32:34.72 


-27:39:58.3 


2002kd 


Prodo 


Dec. 21.64 


03:32:22.34 


-27:44.26.9 


2002ki 


Nanna 


Nov. 22.69 


12:37:28.35 


+62:20:40.0 


2003aj 


Inanna 


Feb. 3.19 


03:32:44.33 


-27:55:06.4 


2003ak 


Gilgamesh 


Feb. 3.19 


03:32:46.90 


-27:54:49.4 


2003az 


Torngasek 


Feb. 20.91 


12:37:19.67 


+62:18:37.5 


2003bd 


Anguta 


Feb. 21.95 


12:37:25.06 


+62:13:17.5 


2003be 


Qiqirn 


Feb. 22.08 


12:36:25.97 


+62:06:55.6 


2003dy 


Borg 


Apr. 4.67 


12:37:09.16 


+62:11:29.0 


20031v 


Vilas 


Apr. 4.67 


12:37:29.00 


+62:11:27.8 


2003eb 


McEnroe 


Apr. 5.65 


12:37:15.18 


+62:13:34.6 


2003eq 


Elvis 


May 24.7 


12:37:48.34 


+62:13:35.3 


2003es 


Ramone 


May 25.5 


12:36:55.39 


+62:13:11.9 
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"standard" rectifications for the camera gain, overscan, spatial bias, dark current, and 
flat fielding. Due to the significant geometric distortion of the ACS WFC (the cost of 
minimizing reflections), we applied the drizzle algorithm (Fruchter & Hook 1997) in the 
Multidrizzle software package (Koekemoer et al. 2004). Because ACS WFC images are 
undersampled at wavelengths shortward of 11,000 A, a better sampled and more precise SN 
PSF can be obtained by "drizzling" (i.e., resampling and combining) the images at a pixel 
scale flner than the physical ACS WFC size of 0.05" pixel""^. However, such improvements 
can only be realized with well-dithered images. The relative size of the dither was measured 
for each frame using source catalogs. Nearly all of the F850LP images in the survey and its 
follow-up were obtained at 4 independent dither positions and were subsequently resampled 
to 0.033" pixel- ^ Imaging in F775W and F606W^ utilized only 2 dither points and the 
physical pixel scale was maintained. 

For NICMOS reductions the CALNICA pipeline in the STSDAS package in IRAF was 
used to provide calibrated frames. Then the wcll-dithered frames were drizzled to half the 
physical pixel scale of Camera 2, i.e., 0.038" pixel"^. The size of each frame's dither was 
determined by cross-correlation of common sources. 

The phasing of GOODS (spanning more than 200 days) provided additional observations 
of the host galaxies with negligible SN light to serve as subtraction templates for the ACS 
data. Subsequent light-curve fltting was used to estimate the expected brightness of the SN 
at the time of each GOODS epoch. Epochs with negligible contamination from SN light 
were combined to obtain deep subtraction templates. By the last epoch of GOODS, SNe 
discovered in the flrst search had faded sufficiently (~2 months past maximum and 3 to 4 
mag below peak) to provide useful templates. For SNe found in all subsequent searches, the 
initial survey epochs provided the components of deep templates. 

Due to the remarkable stability of rectified ACS images it was generally not necessary 
to "blur" (i.e., convolve) images to match the PSFs at different epochs (a conclusion 
reached by fixed-aperture tests on field stars in successive epochs) . Subtractions free from 
host contamination were obtained by fiux-conserving registration of the templates to the 
follow-up frames. The two exceptions were SN 2002hp and SN20031v which resided < 0.05" 
from the sharp nucleus of bright elliptical hosts and for which the technique of matching 
PSFs was used (Alard & Lupton 1998). 

For NICMOS imaging, the necessity of obtaining "clean" subtraction templates was 
judged for each object based on the complexity of the SN site in the bluer, F850LP 
templates. We judged late-time templates were needed and were obtained for SN 2002hp, 



Science Foundation. 
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SN 2002ki, and SN 2003az. As with the ACS data, only SN 2002hp necessitated image 
convolution to match the PSF before subtraction. 

The magnitudes of the SNe in the ACS images were calculated by fitting a PSF 
produced using bright comparison stars scaled to match the infinite-aperture zero-points 
of Sirianni et al. (2004). All ACS passband magnitudes are given as Vega- normalized 
magnitudes. Residual "sky" flux was measured using annuli centered on the SN and with 
an inner and outer radius of 0.66" and 1.00", respectively. The center of the PSF fit was 
determined from a centroid of a stack of all SN images. 

Gilliland & Riess (2002) have shown that the encircled energy of red stars in the IR is 
more dispersed than for blue stars, likely due to backside scattering from the ACS WFC 
CCD mounting. This is an important effect for images in F850LP with SNe la at z > 0.8 
which are considerably red {i — z ^ 1 mag) and can be expected to suffer the "red halo" 
effect. To measure the SNe la in these images, we used a few bright, red {i — z = 1.0 mag) 
field stars as a PSF template. 

Statistical uncertainties for the SN magnitudes were determined by adding and 
recovering artificial PSFs with the measured SN flux (Schmidt et al. 1998). Additional 
uncertainty was included for the shot noise of the measured SN flux. We found excellent 
agreement between the uncertainties calculated empirically and synthetically from the 
STScI exposure-time calculator. 

NICMOS measurements of supernovae were similar, except out model PSF and 
zero-points were calculated from observations of standard stars P330E and G191B2B in 
Cycle 11 (M. Dickinson et al. 2003, private communication). For SN 2003es and SN 2003ak 
which have non-complex backgrounds, the underlying "sky" was determined from isophotal 
modehng of the host. For SN2003dy and SN2003eq the mean background was estimated 
with local apertures. Table 2 includes the measured magnitudes of the SNe la presented 
here. Figure 2 shows their light curves. 

2.3. Spectra 

Supernovae are best classified by the presence and absence of diagnostic features in 
their spectra (see Filippenko 1997 for review). For 14 of the 17 SNe hsted in Table 1, we 
obtained spectra of the SNe near maximum light. For two of these (SN 20031v and SN 
2002hp), bright, elliptical hosts yielded redshifts but overwhelmed the spectra of the SNe. 
For the remaining three (SN 2003ak, SN 2003aj, and SN 2002fw), spectra of the hosts to 
identify their redshifts were obtained when the SNe were no longer visible. Spectra of the 
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twelve visible SNe are shown in Figure 3, and details of all the spectroscopy are provided in 
Table 3. 

To classify the SNe, the detected SN spectra (shown in Figure 3) were cross-correlated 
with template spectra (after removal of the continuum) to identify their type and redshift 
using the "SNID" algorithm (Tonry et al. 2003). For the cases listed in Table 3 for 
which narrow-line host emission was identified, the redshift was constrained to the value 
determined from the host emission before cross-correlation. For cases where the spectra 
were of low S/N, the significance of a SN la classification is much greater if the SN redshift 
is fixed a priori by host emission. For all 12 spectra shown in Figure 3, SNID provided a 
significant classification for each as type la. Although the diagnostic used by the SNID 
algorithm rehes on the whole spectrum, the majority of these SNe can also be classified as 
type la from the presence of Si II absorption at 4130 A (Coil et al. 2001). Specifically, Si II 
absorption is detected in the two highest-redshift spectra presented here, SN 2003dy and 
SN 2002fw. Broad Ca II absorption near 3750 A is visible in all the spectra as well, but this 
feature is much less secure than Si II for SN la classification. 

For two SNe (SN 20031v and SN 2002hp) whose spectra were dominated by those of 
their hosts, the nature of their red, eUiptical hosts allowed us to classify them as highly 
probable SNe la (as was previously the case with SN 1997ff; Riess et al. 2001). 

For the three SNe without any spectroscopy but with host redshifts, classification 
requires greater consideration. Based on the UV color selection method described by Riess 
et al. (2003), SN 2002fx, SN 2003ak, and SN 2003aj are hkely to be SNe la because of 
their red discovery far-UV colors. However, further follow-up of SN 2003aj yielded a rapid 
decline, uncharacteristic of normal SNe la (see Strolgcr et al. 2004). The photometric 
records of the other two arc consistent with SNe la. A pre-discovery observation of SN 
2003ak on the rise is also consistent with the narrow range of SN la rise behavior which 
rise in 3 weeks (Riess et al. 1999b), but not with the majority of core-collapse SNe which 
typically rise in 2 weeks or less, indicating it is very likely to be an SN la. 

For the 16 SNe la (excluding SN 2003aj), we have sufficient quality of photometry to 
yield robust luminosity distances. 
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Table 2 


. SN la Imaging 






Vega Mag 


Epoch(rest) 


K-Corr 




SN 2002fx 




F775W 




F775W U 


495.00 


28.0(1.0) 


-13.5 


0.001(0.06) 


537.80 


27.07(0.25) 


4.2 


0.33(0.05) 


580.00 


28.97(1.00) 


21.8 


0.81(0.02) 




F850LP 




F850LP B 


490.39 


27.48(0.45) 


-15.5 


-1.48(0.02) 


537.79 


25.17(0.07) 


4.2 


-1.16(0.03) 


580.49 


27.07(0.27) 


22.0 


-0.79(0.03) 


SN 2003eq 




F606W 




F606W -> U 


783.68 


24.55(0.01) 


2.7 


0.73(0.03) 




F775W 




F775H^ ^ B 


783.68 


23.19(0.01) 


2.7 


-1.19(0.02) 


799.09 


23.64(0.02) 


11.0 


-1.15(0.02) 


807.27 


24.12(0.04) 


15.4 


-1.12(0.02) 


819.79 


25.02(0.07) 


22.2 


-0.99(0.02) 


838.30 


26.18(0.30) 


32.2 


-0.96(0.02) 




F850LP 




F850LP V 


783.68 


23.03(0.01) 


2.7 


-1.28(0.02) 


792.10 


23.18(0.02) 


7.2 


-1.22(0.02) 


799.09 


23.37(0.02) 


11.0 


-1.17(0.03) 


807.27 


23.72(0.04) 


15.4 


-1.12(0.06) 


819.79 


24.19(0.06) 


22.2 


-0.93(0.05) 


838.30 


24.95(0.07) 


32.2 


-0.69(0.02) 




FllOW 




FllOW 


792.94 


23.36(0.10) 


7.7 


-1.24(0.02) 


SN 2003es 




F775W 




F775W U 


784.50 


24.09(0.03) 


7.7 


-0.87(0.02) 




i^850LP 




F850LP B 


784.50 


23.67(0.02) 


7.7 


-1.35(0.02) 


792.37 


24.05(0.05) 


11.7 


-1.37(0.02) 


801.27 


24.64(0.06) 


16.2 


-1.43(0.03) 


807.81 


24.90(0.07) 


19.6 


-1.49(0.03) 


820.91 


25.76(0.08) 


26.2 


-1.55(0.02) 


838.10 


26.07(0.15) 


35.0 


-1.56(0.02) 




FllOW 




FllOW ^ V 


792.76 


24.20(0.08) 


11.9 


-1.32(0.02) 


SN 2003az 




F775W 




F775W U 
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Table 2 — Continued 





Vega Mag 


Epoch(rest) 


K-Corr 


690.89 


25.06(0.05) 


0.0 


-0.22(0.03) 


701.16 


25.42(0.05) 


4.5 


-0.15(0.02) 




F850LP 




F850LP B 


690.89 


24.31(0.04) 


0.0 


-1.43(0.03) 


701.16 


24.44(0.04) 


4.5 


-1.39(0.02) 


709.09 


24.60(0.05) 


8.0 


-1.38(0.02) 


716.92 


25.04(0.06) 


11.5 


-1.37(0.02) 


726.52 


25.47(0.08) 


15.7 


-1.33(0.03) 


733.25 


25.71(0.09) 


18.7 


-1.27(0.04) 




FllOW 




FllOW 


703.62 


24.22(0.06) 


5.6 


-1.59(0.03) 


710.58 


24.37(0.06) 


8.7 


-1.54(0.03) 


SN 2002kc 




F606W 




F&Q&W V 


629.62 


22.27(0.01) 


-7.5 


-0.23(0.02) 


672.33 


23.25(0.01) 


27.5 


0.23(0.02) 




F775W 




F775W R 


629.62 


21.78(0.02) 


-7.5 


-0.43(0.02) 


642.50 


21.32(0.10) 


3.0 


-0.47(0.02) 


672.33 


22.08(0.01) 


27.5 


-0.44(0.02) 




F850LP 




FS50LP I 


629.62 


21.66(0.01) 


-7.5 


-0.46(0.02) 


672.33 


21.89(0.01) 


27.5 


-0.19(0.04) 


SN 2003eb 




F606W 




F606W U 


734.58 


24.23(0.02) 


-1.3 


0.84(0.04) 


783.47 


27.02(0.15) 


24.0 


1.36(0.02) 




F775W 




F77bW B 


734.65 


23.02(0.02) 


-1.3 


-1.15(0.02) 


745.65 


23.12(0.02) 


4.3 


-1.11(0.02) 


783.54 


25.18(0.06) 


24.1 


-0.83(0.02) 


799.10 


25.87(0.10) 


32.2 


-0.79(0.02) 




F850LP 




F850LP ^ V 


734.65 


22.79(0.02) 


-1.3 


-1.34(0.02) 


745.65 


22.81(0.02) 


4.3 


-1.29(0.05) 


751.12 


22.94(0.02) 


7.2 


-1.22(0.04) 


763.61 


23.45(0.02) 


13.7 


-1.11(0.08) 


773.89 


23.98(0.04) 


19.1 


-0.91(0.08) 


783.54 


24.42(0.04) 


24.1 


-0.72(0.06) 


792.37 


24.72(0.05) 


28.7 


-0.59(0.02) 


792.10 


24.79(0.06) 


28.5 


-0.59(0.02) 
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JD" 


Vega Mag 


Epoch(rest) 


K-Corr 


799.09 


25.05(0.08) 


32.2 


-0.52(0.02) 


801.28 


24.93(0.06) 


33.3 


-0.52(0.02) 


807.81 


25.07(0.07) 


36.7 


-0.56(0.02) 


820.92 


25.33(0.07) 


43.5 


-0.60(0.02) 


838.10 


25.37(0.08) 


52.5 


-0.66(0.02) 


SN 20031v 




F606W 




F606W U 


733.65 


25.14(0.05) 


6.0 


0.857(0.02) 




F775W 




i^775M^ ^ B 


733.65 


23.54(0.05) 


6.0 


-1.12(0.02) 


745.65 


24.09(0.05) 


12.5 


-1.09(0.02) 


783.54 


26.36(0.25) 


32.8 


-0.91(0.02) 




F850LP 




i^850LP ^ V 


692.42 


26.39(0.30) 


-16.1 


-1.17(0.02) 


733.65 


23.25(0.04) 


6.0 


-1.23(0.03) 


745.65 


23.57(0.05) 


12.5 


-1.13(0.04) 


751.12 


23.86(0.06) 


15.4 


-1.05(0.07) 


763.61 


24.60(0.08) 


22.1 


-0.80(0.05) 


773.89 


24.92(0.12) 


27.6 


-0.68(0.02) 


783.54 


25.63(0.15) 


32.8 


-0.65(0.02) 


799.09 


26.05(0.25) 


41.2 


-0.71(0.02) 


807.27 


26.28(0.30) 


45.6 


-0.73(0.02) 


SN 2002hr 




i^606W^ 




F606H^ ^ U 


579.64 


24.01(0.03) 


-6.3 


0.09(0.04) 


629.47 


25.78(0.15) 


26.2 


-0.38(0.02) 


674.11 


27.38(0.30) 


55.5 


-0.37(0.04) 




F775W^ 




F77bW B 


579.64 


23.53(0.03) 


-6.3 


-0.86(0.02) 


589.68 


23.17(0.12) 


0.2 


-0.89(0.02) 


590.01 


23.25(0.11) 


0.4 


-0.89(0.02) 


596.65 


23.41(0.14) 


4.7 


-0.95(0.05) 


614.50 


23.88(0.08) 


16.4 


-1.29(0.11) 


629.47 


24.39(0.05) 


26.2 


-1.68(0.02) 


674.11 


25.62(0.10) 


55.5 


-1.60(0.02) 




F850LP 




F850i:P ^ V 


579.64 


23.34(0.03) 


-6.3 


-0.97(0.02) 


629.47 


23.75(0.04) 


26.2 


-1.19(0.02) 


674.11 


24.97(0.07) 


55.5 


-1.15(0.02) 



SN 2003bd 



Fmm U 
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Vega Mag 


Epoch(rest) 


K-Corr 


691.94 


24.54(0.04) 


11.6 


0.161(0.02) 


735.44 


27.45(0.12) 


37.6 


0.180(0.02) 








Fll^W B 


691.94 


23.38(0.03) 


11.6 


-1.05(0.02) 


735.44 


25.72(0.10) 


37.6 


-1.26(0.02) 


745.65 


26.00(0.10) 


43.7 


-1.25(0.02) 




F850LP 




F850LP V 


642.50 


26.98(0.40) 


-17.9 


-1.10(0.02) 


691.94 


23.14(0.03) 


11.6 


-1.08(0.02) 


735.44 


24.75(0.05) 


37.6 


-1.03(0.02) 


745.65 


24.95(0.06) 


43.7 


-1.03(0.02) 


751.18 


25.04(0.07) 


47.0 


-1.03(0.02) 


763.60 


25.30(0.08) 


54.5 


-1.04(0.02) 


773.95 


25.40(0.08) 


60.7 


-1.04(0.02) 


792.09 


25.72(0.10) 


71.5 


-1.05(0.02) 


807.27 


25.67(0.25) 


80.6 


-1.05(0.02) 


819.79 


25.72(0.25) 


88.1 


-1.05(0.02) 


SN 2002kd 








F606W -> U 


629.42 


24.93(0.03) 


-8.9 


0.326(0.02) 


673.52 


25.78(0.05) 


16.4 


0.330(0.02) 








F775T^ ^ B 


629.43 


23.75(0.02) 


-8.9 


-1.07(0.02) 


644.50 


22.93(0.17) 


-0.2 


-1.11(0.02) 


645.50 


22.73(0.10) 


0.3 


-1.11(0.02) 


673.53 


24.15(0.05) 


16.4 


-1.12(0.02) 








F850LP V 


629.47 


23.75(0.01) 


-8.9 


-1.11(0.02) 


639.38 


23.04(0.02) 


-3.1 


-1.13(0.02) 


673.58 


23.73(0.01) 


16.5 


-1.01(0.02) 


SN 2003be 




Fmm 




F606W U 


692.04 


24.44(0.03) 


12.1 


0.087(0.02) 


732.46 


27.21(0.20) 


36.8 


0.091(0.02) 




F775T^ 




F775W B 


641.30 


28.0(0.8) 


-18.7 


-0.72(0.06) 


692.04 


23.42(0.03) 


12.1 


-1.05(0.03) 


732.47 


25.46(0.10) 


36.8 


-1.29(0.02) 


784.28 


26.04(0.10) 


68.4 


-1.21(0.02) 




ii'850LP 




F850LP V 


641.30 


28.0(0.8) 


-18.7 


-1.05(0.02) 
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Vega Mag 




Epoch(rest) 


K-Corr 


692.04 


23.04(0.02) 




12.1 


-1.07(0.02) 


732.46 


24.52(0.03) 




36.8 


-1.09(0.02) 


784.28 


25.53(0.10) 




68.4 


-1.07(0.02) 


SN 2003dy 




F850LP 






F850LP U 


692.42 


27.47(0.75) 




-17.1 


-0.86(0.04) 


733.65 


24.43(0.05) 




0.3 


-0.96(0.04) 


745.65 


24.62(0.07) 




5.4 


-1.05(0.02) 


751.11 


24.75(0.07) 




7.7 


-1.06(0.02) 


763.60 


25.22(0.08) 




12.9 


-1.12(0.02) 


773.89 


25.53(0.10) 




17.3 


-1.16(0.02) 


783.54 


26.87(0.34) 




21.4 


-1.17(0.02) 


801.28 


27.00(0.43) 




28.9 


-1.21(0.02) 


807.81 


27.21(0.64) 




31.6 


-1.19(0.02) 




FllOW 






FllOW B 


751.59 


24.50(0.07) 




7.9 


-1.71(0.02) 


754.59 


24.60(0.08) 




9.1 


-1.73(0.02) 




F160W^ 






Fimw -> R 


751.72 


23.90(0.08) 




7.9 


-1.95(0.02) 


SN 2002ki 




F775W 






F775W ^ U 


600.76 


24.76(0.10) 




-0.3 


-0.47(0.02) 


643.57 


26.78(0.17) 




19.6 


-0.48(0.02) 




F850LP 






F850LP B 


600.76 


23.85(0.04) 




-0.3 


-1.50(0.02) 


643.57 


25.78(0.10) 




19.6 


-1.38(0.02) 


652.21 


25.98(0.10) 




23.7 


-1.34(0.02) 


663.68 


26.88(0.20) 




29.0 


-1.34(0.02) 




F160W 






Fimw I 


652.44 


24.77(0.25) 




23.8 


-1.65(0.02) 


664.47 


24.96(0.25) 




29.4 


-1.55(0.03) 


SN 2003ak 




F850LP 






F850LP ^ U 


627.35 


27.33(0.40) 




-14.0 


-0.72(0.02) 


673.13 


25.63(0.08) 




3.8 


-0.70(0.02) 


680.12 


25.82(0.08) 




6.6 


-0.69(0.02) 


694.48 


26.52(0.12) 




12.2 


-0.65(0.02) 




FllOW 






FllOW ^ B 


681.06 


25.15(0.10) 




7.0 


-1.72(0.02) 


693.01 


25.25(0.10) 




11.6 


-1.75(0.02) 


708.48 


25.60(0.15) 




17.7 


-1.82(0.04) 
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JD" 


Vega Mag 




Epoch(rest) 


K-Corr 


715.35 


26.02(0.15) 




20.4 


-1.87(0.02) 




Fmow 






F160W V 


681.26 


24.10(0.05) 




7.0 


-2.32(0.02) 


693.41 


24.30(0.07) 




11.8 


-2.28(0.02) 


701.61 


24.70(0.08) 




15.0 


-2.25(0.07) 


SN 2002hp 




F775W 






F775W U 


579.49 


25.49(0.10) 




3.2 


-0.05(0.04) 




F850LP 






FSbOLP B 


537.31 


26.80(0.39) 




-15.0 


-1.50(0.03) 


579.49 


24.29(0.04) 




3.2 


-1.30(0.02) 


589.05 


24.86(0.06) 




7.3 


-1.26(0.02) 


595.38 


25.04(0.08) 




10.1 


-1.25(0.02) 


603.81 


25.65(0.08) 




13.7 


-1.22(0.02) 


613.80 


26.07(0.14) 




18.1 


-1.13(0.05) 


629.30 


26.59(0.17) 




24.8 


-1.02(0.02) 


639.38 


27.29(0.35) 




29.2 


-1.01(0.02) 




FllOW 






FllOW V 


589.17 


24.29(0.07) 




7.4 


-1.54(0.03) 


595.50 


24.44(0.07) 




10.1 


-1.49(0.03) 


SN 2002fw 




F775W 






F77hW U 


536.80 


25.34(0.05) 




-5.8 


-0.24(0.02) 


548.30 


24.98(0.08) 




-0.8 


-0.16(0.04) 


578.40 


26.41(0.10) 




12.2 


0.09(0.04) 




i^850i:P 






F850LP B 


536.80 


24.57(0.05) 




-5.8 


-1.43(0.02) 


548.30 


24.25(0.06) 




-0.8 


-1.43(0.03) 


552.30 


24.47(0.06) 




0.8 


-1.43(0.03) 


557.40 


24.44(0.06) 




3.0 


-1.39(0.02) 


567.80 


24.56(0.08) 




7.6 


-1.36(0.02) 


577.50 


24.99(0.06) 




11.8 


-1.34(0.02) 


578.50 


25.13(0.07) 




12.2 


-1.33(0.02) 


595.40 


25.97(0.16) 




19.6 


-1.19(0.04) 


603.40 


26.88(0.25) 




23.0 


-1.09(0.02) 


628.00 


27.31(0.30) 




33.7 


-1.06(0.02) 




FUOW 






FllOW ^ V 


549.51 


24.15(0.08) 




-0.3 


-1.69(0.02) 


557.57 


24.25(0.08) 




3.1 


-1.66(0.03) 




F160W 






F160W ^ R 


549.75 


23.92(0.08) 




-0.2 


-1.92(0.05) 
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3. Light-Curve Fitting 

Distances to SNe la with individual precision approaching 7% can be measured utihzing 
empirical relationships between light-curve shape and luminosity as well as color-curve 
shape and extinction. Our primary analysis uses a revision of the multi-color light-curve 
shape (MLCS) fitting method (Riess, Press, & Kirshner 1995, 1996a; Riess et al. 1998) 
described by Jha (2002) and Jha et al. (2004a), and hereafter referred to as "MLCS2k2." 
Previous versions of MLCS fit rest-frame BVRI, but MLCS2k2 includes U band templates 
based on a new set of 25 well-observed SNe la in the U band from Jha (2002) and Jha 
et al. (2004b). This UV extension of MLCS allows us to extend the Hubble diagram of 
SNe la to 2; > 1 using light curves observed in the reddest available band on ACS. The 
results from Jha (2002) demonstrate that rest-frame f/-band light curves provide similar 
information on the epoch of maximum, relative luminosity, reddening, and distance as 
optical hght curves, albeit with lower precision (by a factor of 1.5). Because the MLCS2k2 
covariance matrix contains the bandpass-dependent variance of SNe la (primarily in the 
form of autocorrelation along the diagonals as well as in the form of two-point covariance in 
off-diagonal terms), the [/-band light curves can be used together with optical light-curve 
data with the appropriate weight and propagation of uncertainty in the multi-color fit. 

Additional improvements to MLCS2k2 include a more self-consistent treatment of the 
allowed range of extinction and extinction laws, as well as an improved determination of 
the unreddened SN la color (see Jha et al. 2004a for details). This method still empirically 
models a light curve as the sum of a fiducial template and a set of phase-dependent vectors 
(linear and quadratic) whose contribution scales with the luminosity offset from the fiducial 
curve at peak. Luminosity corrections are not extrapolated beyond the range observed 
in the local sample. Simultaneous fitting in multiple colors constrains the line-of-sight 
reddening and distance. 



Table 2 — Continued 







Vega Mag 


Epocli(rest) 


K-Corr 


557.76 




23.85(0.09) 


3.2 


-1.86(0.03) 


Actually JD- 
Uncertainties in 


-2,450,000. 
magnitudes 


are listed in parentheses. 







Table 3. Spectroscopic Data 



SN 


UT Dat(- 


Iiistmiiicnt 


exposiirefsec) 







2002fw 


Sep. 31, 2002 


HST ACS 


15000 


1.30'^'! 


2002fx 


Scp.14, 2003 


Keck-II NIRSPEC 


2000 




2002hp 


Nov. 7, 2002 


Keck-I LRIS 


7800 


1.305'''2 




Nov. 7, 2002 


VLT FORS 


14000 


1.305'''2 


2002hr 


Nov. 8, 2002 


Kcck-I LRIS 


7800 


0.526'='! 


2002kc 


Jan. 7, 2003 


Kcck-I LRIS 


1500 


0.216'='! 


2002kd 


Jan. 1, 2003 


Magellan LDSS 


7200 


0.735'='! 


2002ki 


Jan. 7, 2003 


Keck-I LRIS 


2700 


1.141'='! 


2003aj 


Oct 1-3, 2003 


VLT F0RS2 


16800 


1.307'''^ 


2003ak 


Sep. 11, 2003 


Keck-II NIRSPEC, VLT F0RS2 


14000 


1.551'''3 


2003az 


Mar. 3, 2003 


HST ACS 


6500 


1.27°'! 


2003bd 


Feb. 27/28, 2003 


Keck-I LRIS 


16500 


0.67°'! 


2003bc 


Feb. 28, 2003 


Keck-I LRIS 


5400 


0.64'='! 


2003dy 


Apr. 16, 2003 


HST ACS 


15000 


1.34'*'! 


20031v 


Apr. 16, 2003 


HSTkCS 


15000 


0.935'^'2 


2003eb 


Apr. 16, 2003 


HST ACS 


15000 


O.SQQ"^'! 


2003eq 


Jun. 2, 2003 


HST ACS 


6000 


0.839'^'! 


2003es 


Jun. 2, 2003 


HST ACS 


6000 


0.954'''! 



^From cross-correlation with broad SN features. 
!'From narrow features in the host-galaxy spectrum. 
'=Froni both (a) and (b). 

•^From (a) and Cowie et al. (2004), Wirth et al (2004). 
! Classified as SN la with high confidence from spectrum. 
^Classified as SN la with high confidence from early-type, red host. 
^Photometric properties indicate likely SN la. 
^Uncertain type. 

For VLT host spectra of 2003ak,aj see Nonino et al 2004 
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Extinction priors are used together with the observed reddening to estimate the 
expected dimming of the SN magnitudes resulting from dust. Recent analyses by the high-z 
supernova search team (HZT; Schmidt et al. 1998) and the supernova cosmology project 
(SCP; Perlmutter et al. 1997), published respectively by Tonry et al. (2003) and Knop et 
al. (2003), make use of a Galactic prior for the extinction law (which is corroborated at 
low redshifts; Riess, Press, Kirshner 1996b) and flag SNe with large reddening to reduce 
systematic errors from non-Galactic-type dust. The HZT has employed an exponential 
extinction prior whose functional form derives from modeled lines-of-sight (Hatano, Branch, 
& Deaton 1998) and from the a posteriori distribution of extinction values. Here we utilize 
the exponential prior on extinction and the gaussian prior on the intrinsic color derived 
from the a posteriori distribution of Jha et al. (2004a) and flag SNe with large measured 
reddening as unreliable. Each SN la is corrected for galactic reddening as estimated in the 
direction of each SN by Schlegel, Finkbeiner, & Davis (1999) before the observed colors are 
used to estimate the host reddening. 

K-corrections are used to account for the SN redshift and provide a transformation 
between an observed-frame magnitude and a rest-frame magnitude (Okc & Sandage 1968). 
We use composite spectra of SNe la from Nugent et al. (2003) to calculate Vega-normalized 
"cross-band" K-corrections. Individual K-corrections were calculated to the best-matching 
passbands for the appropriate phase and redshift of the SN. Each SN is then fit using 
MLCS2k2 to provide a custom model of the phase-dependent colors. The colors of the 
Nugent et al. (2003) spectral energy distributions are then matched to the model by 
multiplication of a spline-interpolation. Next, the K-corrections are recalculated, and 
the process continues until convergence (usually after two to three iterations). The final 
K-corrcctions are given for each measured magnitude in Table 2. The MLCS2k2 fits to each 
of the 16 new SNe la are shown in Figure 2 and the fit parameters are given in Table 4. 

We have also used an additional hght-curve fitting method, "Bayesian Template 

Method" (BATM; Tonry et al. 2003), to provide independent estimates of the luminosity 
distances of the new SN la data presented here. BATM is a "template-fitting" method 
which seeks to identify close matches between an individual SN la and a well-observed, 
local counterpart and then calculates their distance ratio by synthetically redshifting the 
nearby template to the observed frame. The current realization of BATM has far fewer 
[/-band examples than MLCS2k2, because BATM was developed without the Jha (2002) 
data. Consequently, BATM did not converge on a solution for three sparsely observed 
SNe la (SNe 2002kc, 2003be, and 2002fx), whose MLCS2k2 fits relied on sampling in the 
rest-frame U band. 
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3.1. The Ground-Based Discovery Set 

We construct the expansion history of the Universe by using this new set of HST- 
discovered objects together with pubUshed observations of supernovae over a wide range 
in redshift. Tonry et al. (2003) have recently compiled the distances and redshifts for 172 
SNe la. Using available results from different light-curve fitting methods (including BATM, 
MLCS, Amis, snapshot, and stretch) for each SN la, Tonry et al. (2003) corrected for 
zeropoint differences between methods and provided best estimates of the distance to each 
SN la from a median of the distance estimates from individual methods. 

Although the Tonry et al (2003) data set represented the state of the art in February, 
2003, when it was submitted, there have been some significant developments since then that 
need to be included to build the most relaible data set for analyzing the HST-discovered 
objects. For the SCP, Knop et al. (2003) report on a new set of 11 SNe la at OA < z < 0.85 
as well as on a reanalysis of the original high-redshift SNe la from the SCP (Perlmutter 
et al. 1999). In this reanalysis they now exclude 15 of the 42 high-redshift SNe from 
Perlmutter et al. (1999) due to inaccurate color measurements and uncertain classification. 
Knop et al. (2003) flag an additional 6 of the original 42 SNe, as well as 5 of the new 
11 SNe, as likely SNe la but faihng a "strict SN la" sample cut. For the HZT, Harris et 
al. (2004) report on a large set of new high-redshift SNe (22 in all), with widely varying 
degrees of completeness of the spectroscopic and photometric records. Finally, Blakeslee et 
al. (2003) report on 2 new SNe la discovered with ACS on HST. The development of the 
MLCS2k2 method, which includes [/-band observations when they are available, makes it 
worthwhile to revisit the previously published data. 

We recompiled a set of previously observed SNe la relying on large, published samples, 
whenever possible, to reduce systematic errors from differences in calibration. To compile 
SNe la at 0.01 < z < 0.15 we used the 3 largest, modern data sets of such SNe la published 
to date: the Calan-Tololo Survey (29 SNe la; Hamuy et al. 1996), the CfA Survey I (22 
SNe la; Riess et al. 1999a), and the CfA Survey II (44 SNe la; Jha et al. 2004b). ^ At 
higher redshifts we used SN la photometry published by Riess et al. (1998, 2001), Tonry 
et al. (2003), Knop et al. (2003), and Barris et al. (2004), as well as SNe la tabulated by 
Tonry et al. (2003) from Suntzeff et al. (2004), Leibundgut et al. (2004), Clochiatti et al. 
(2004), and Jha et al. (2004c). Despite the apparently large number of sources of data, 
the majority of the data at z > 0.1 comes from the HZT (Schmidt et al. 1998; Riess et 
al. 1998), and thus the methods used to calibrate all of these SNe la are extremely similar 



^Note that the CfA SNe were generally discovered by other searches, especially the Lick Observatory SN 
Search with the Katzman Automatic Imaging Telescope (Filippenko et al. 2001; Filippenko 2003). 
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(and familiar to the authors of the present paper). 

In order to reduce systematic errors which arise from differences in hght-curve fitting 
(and K-correcting) methods, we have made every effort to consistently refit all the past 
data with a single method, MLCS2k2 (Jha et al. 2004a). An exception was made for the 
high-redshift SNe la from the SCP (Perlmutter et al. 1999) because photometry of these 
SNc la remain unpublished. For these we utilized the reanalyzed distances as given by 
Knop ct al. (2003). We transformed them to the MLCS2k2 distance scale by solving for the 
SN la luminosity zero-point required to match the mean distances to low-redshift objects 
fit by the two methods. 

Unfortunately, a large variation exists in the quality and breadth of the photometric 
and spectroscopic records of individual SNe. Ideally, each SN la would have the same 
well-defined spectroscopic features (Filippenko 1997), a spectroscopic redshift, and 
well-sampled light curves and color curves. However, this is often not the case. As with the 
photometry, most of the extant high-redshift spectra from the SCP have not been published 
(e.g., Perlmutter et al. 1999; Knop et al. 2003). This makes it difficult to apply a uniform 
set of criteria to SNe la in constructing a cosmological sample. 

To reflect the differences in the quality of the spectroscopic and photometric record for 
individual supernovae, we divide the objects into "high-confidence" (hereafter "gold" ) and 
"likely but not certain" (hereafter "silver" ) subsets. Ideally, we would assign each supernova 
a weight in any overall fit that refiected its individual uncertainty. However, distance 
errors resulting from spurious problems more common to lower-confidence SNe la such as 
SN misclassification, large extinction (amplifying uncertain extragalactic extinction laws), 
and poorly constrained colors are difficult to quantify. So we use the coarser approach of 
separating the high-confidence gold events from the larger set. 

We adopted in our gold set all SNe la which were included by the original sources in 
their most stringent subsets (when such discriminations were made). Any SN la fiagged as 
having a cause for a specific concern was not included in this set. The two primary reasons 
for rejecting a SN la from this set are that (1) the classification, though plausible, was not 
compelling (see discussion), and (2) the photometric record is too incomplete to yield a 
robust distance (i.e., the number of model parameters is roughly equal to the number of 
effective samphngs of the light curve) . 

SNe la included in previous cosmological samples but rejected from our gold sample 
include SN 1999fh (poorly constrained light curve; Tonry et al. 2003), SN 1997ck (poor 
color information; Garnavich et al. 1998), all "snapshot" SNe la from Riess et al. (1998), 
15 SNe la from Perlmutter et al. (1999) later discarded by Knop et al. (2003) as well as 11 
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additional SNe flagged by Knop et al. (2003; flag values 1, 2, 3), SNe la from Barris et al. 
(2003) without SNID classification, and any SN la with more than 1 mag of extinction or 
whose light curve begins more than 10 rest-frame days after maximum as determined from 
the MLCS2k2 fit. 

The same criteria were apphed to the GOODS SNe la whose individual classifications 
were described in §2.3. As a result, two of these SNe were rejected from the gold sample: 
SN 2002fx, whose classification was not certain enough, and SN 2002kc, whose fit indicated 
> 1 mag extinction. The gold set contains a total of 157 SNe la. 

The silver set contains the objects identified by the above sources as likely SNe la but 
failing one criterion for inclusion in the gold category. The silver set includes a total of 29 
SNc la. SNe failing more than one criterion were excluded from the analyses. The final 
membership rosters of the subsets are tabulated in the Appendix. 

For most of our cosmological analyses wc focus on results derived from the gold set, but 
for a few analyses with the largest number of free parameters (and thus the most limited 
in statistical inference) we include the silver set (with the caveats arising from its reduced 
reliability) . 



Distance estimates from SN la light curves are derived from the luminosity distance. 



where £ and T arc the intrinsic luminosity and observed flux of the SN within a given 
passband, respectively. Equivalently, logarithmic measures of the flux (apparent magnitude, 
m) and luminosity (absolute magnitude, M) were used to derive extinction-corrected 
distance moduli, /Iq — m — M — Slogd^ -|- 25 (d^ in units of megaparsecs) . In this context, 
the luminosity is a "nuisance parameter" whose value is unimportant for kinematic (and 
most cosmological) studies. We have used the MLCS2k2 method and the data described in 
§2 to derive accurate and individual relative distance moduli for the sets of SNe described 



In Figure 4 we show the Hubble diagram of distance moduh and redshifts for the new 
HST-discovered SNe la in the gold and silver sets. Although these new SNe la span a 
wide range of redshift (0.21 < z < 1.55), their most valuable contribution to the SN la 
Hubble diagram is in the highest-redshift region where they effectively delineate the range 



4. Cosmological Constraints 




(1) 



in §3. 
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at 0.85 < z < 1.55 with 11 new SNe la, including 6 of the 7 highest-redshift SNe known 
(the seventh being SN 1997ff; Riess et al. 2001). 

The relationship between distance and redshift over a significant fraction of the Hubble 
time can be considered either empirically as a record of the (integrated) expansion history 
of the Universe, or theoretically as constraints on the mass-energy terms contained in the 
Friedman equation and affecting the expansion. In the next subsections we consider both 
approaches. 



4.1. Expansion History: A Kinematic Description 

It is valuable to consider the distance-redshift relation of SNe la as a purely kinematic 
record of the expansion history of the Universe, without regard to its cause. An empirical 
description of the time variation of the scale factor, a{t), can provide answers to basic 
questions (e.g., "When was the Universe (if ever) accelerating or decelerating?") and 
model-independent constraints with which to test cosmological models. 

Following Turner & Riess (2002), we empirically define the luminosity distance in 
Euclidean space (i.e., f^totai — 1-0; as motivated by infiation) as the integral of the inverse 
of the preceding expansion rate. 



d, 



c(l + z) 



/■^ du 
Jo H{u) 



^cil + z 



Jo 



exp 



- r [1 + q(u)]dln(l + u)] du, (2) 

^0 J 



where 



H{z) = -, 
a 



(3) 



q{z) ^ {-d/a)/H\z) = ^-^^ - 1. (4) 

Note that equation (2) is not an approximation but is an exact expression for the luminosity 
distance in a geometrically fiat Universe (though generalizable for non-zero curvature), 
given an expression for the epoch-dependent deceleration parameter, q{z), and the present 
Hubble constant, Hq. Here we employ equation (2) as a kinematic model of the SN la data 
with parametric representations for q{z). 
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Given evidence that the Universe has recently been accelerating [i.e., q{z ~ 0) < 0], 
hints that it may have once been decelerating [i.e., q{z > 1) > 0; Riess et al. 2001; Turner 
& Riess 2002], and the large leverage in redshift of the current SN sample, we consider 
resolving q[z) into two distinct components or epochs. A linear two-parameter expansion 
for q{z) which is continuous and smooth is q{z) = qo + zdq/dz, where dq/dz is defined to be 
evaluated at z = 0. 

The likelihood for the parameters q^ and dq/dz can be determined from a statistic, 
where 



, cr^ is the dispersion in supernova redshift (transformed to units of distance moduli) due to 
peculiar velocities and cr^g . is the uncertainty in the individual distance moduli. This term 
also includes the uncertainty in galaxy redshift. Due to the extreme redshift of our distant 
sample and the abundance of objects in the nearby sample, our analysis is insensitive to 
the value we assume for cr„ within its likely range of 200 km s~^ < cr^ < 500 km s~^. For 
our analysis we adopt = 400 km s~^. For high-rcdshift SNc la whose redshifts were 
determined from the broad features in the SN spectrum, we add 2500 km s~^ in quadrature 
to ay. 

Marginalizing our likelihood functions over the nuisance parameter, Hq (by integrating 
the probability density P oc e~'^^/^ for all values of Hq), yields the confidence intervals 
shown in Figure 5. As shown, both the gold set or the gold and silver sets strongly favor 
a Universe with recent acceleration (go < 0) and previous deceleration {dq/dz > 0) with 
99.2% and 99.8% likelihood (summed within this quadrant), respectively. With this same 
model we can also derive the likelihood function for the transition redshift, Zt, defined as 
q{zt) — 0. Summing the probability density in the qo vs. dq/dz plane along lines of constant 
transition redshift, z^ — —qo/ (dq/dz), yields the likelihood function in Figure 5. We find a 
transition redshift of Zf = 0.46 ± 0.13. In Figure 6 we show the Hubble diagram for the 
SNe la compared to a discrete set of kinematic models. 

An alternate, kinematic model is derived using the first three time derivatives of the 
scale factor. Following Visser (2003), the Hubble, deceleration, and jerk parameters are 
defined as 



X^{Ho,qo, dq/dz) = ^ 




(5) 



H{t) = +d/a , 



(6) 



{d/a){d/a) ^ , and 



(7) 
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j(t) = +(a/a)(d/a)-^ (8) 

The deceleration and jerk parameters are dimensionless, and a Taylor expansion of the 
scale factor around to provides 



a{t) = ao {l + Ho {t - to) - ^ qo H^ {t - to)^ + ^ Jo Hi {t - tof + 0{[t - to]')} , 

(9) 

and hence for the luminosity distance (in Euclidean space), 

dL{z) = ^ {l + ^ [1 - go] ^ - ^ [l - go - 3%' + Jo] + 0{z')^ (10) 
(of. Visser 2003). 

Though related, the jo parameter as defined here and by Visser (2003) is not precisely 
equivalent to our previous dq/dz parameter, providing an alternative parameterization. 
The SN subsets constrain the jo parameter to the positive domain at the 92% to the 95% 
confidence level. That is, the expansion history over the range of the SN data is equally 
well-described by recent acceleration and a constant jerk. Models with discrete values of jo 
are shown in Figure 6. 

In summary, we find strong evidence for a change in the sign of cosmic acceleration in 
the past. 



4.2. Cosmological Constant or Astrophysical Dimming? 

SNe la at 2; ~ 0.5 appear fainter by ~0.25 mag relative to a Universe with VLm = 0.3 
and = 0, a result readily accommodated by a cosmological constant with JIa = 0.7 
(Riess et al. 1998; Perlmutter et al. 1999). Despite the lack of any independent evidence, 
an alternative explanation for this dimming could lie in the astrophysics of supernovae or 
in the propogation of their light to us. Speculative models for astrophysical contamination 
of the SN la signal have been posited; these include extragalactic gray dust with negligible 
tell-tale reddening or added dispersion (Aguirre 1999a,b; Rana 1979, 1980), and a pure 
luminosity evolution (Drell, Lorcdo, & Wasscrman 2000). Here we limit our consideration 
to the observable differences between these hypotheses. 
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The luminosity distance expected in a Priedmann-Robertson- Walker (FRW) cosmology 
with mass density flM and vacuum energy density (i.e., the cosmological constant) JIa is 



dL = cH^\l + z) \nk\~'^^ smn{\nk\'^^ / dz[{l + zf{l + ^mz) - z{2 + z)nA\-^/^], (11) 



where Q,k — 1 — — ^a, and "sinn" is sinh for Qj. > and sin for Qfe < (Carroll, Press, 
& Turner 1992). For = 0, equation (11) reduces to cHq^{1 -\- z) times the integral. With 
dL in units of megaparsecs, the predicted distance modulus is 



Following Goobar, Bergstrom, & Mortsell (2002) we consider two models of gray 
extinction by a homogeneous component of dust: Pdnst{z) = Pdust(l + ^)") where 



The "high-z dust" model represents a smooth background of dust present (presumably 
ejected from galaxies) at a redshift which is greater than the SN sample (i.e., z > 2) and 
diluting as the Universe expands. The total extinction is then calculated as the attenuation 
integrated along the photon path, —2.5log{e:K.p{jQPdust{z)r{z)dz)), where r{z) is the 
coordinate distance traversed by the SN photons. A single free (opacity) parameter is fixed 
by requiring the total extinction at z ~ 0.5 to match the observed peak brightness of SNe la 
in a cosmology with fl\ — 0. 

The "replenishing dust" represents a constant density of dust which is continually 
replenished at precisely the same rate in which it is diluted by the expanding Universe (i.e., 
a = 0). This model is also tuned to match the extinction implied by SNe la at 2; ~ 0.5 in 
the absence of a cosmological constant, hence it requires the tuning of two parameters (as 
well as fast-moving dust which quickly provides a homogeneous background without added 
dispersion from uneven lines-of-sight). We also consider a third model (following Filippenko 
& Riess 2001, Riess et al. 2001, and Blakeslee et al. 2003) to mimic simple evolution which 
scales as z in percent dimming. Our set of models for astrophysical dimming is not an 
exhaustive set of all possibilities, but rather is drawn from physically motivated hypotheses 
(in contrast to purely parametric models of astrophysical dimming, e.g., Drell et al. 2000). 

In Figure 7 we show the Hubble diagram of SNe la relative to the cosmological 
and astrophysical hypotheses. As seen in Table 4, the SN dataset is consistent with an 
Qm — 0.27, Qa = 0.73 cosmology, yielding — 178 for 157 SNe la (degrees of freedom, dof; 
Xdof — 1-13) in the gold set. The total is significantly worse for the high- redshift gray 



Pp = 5\ogdL + 25. 



(12) 




3 for all z "high — z dust," 

for ^ > 0.5 (3 for lower z). "replenishing dust. 
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dust model (Ax^ = 122; 11a for 1 dof) as well as for the simple model of evolution with 
dimming oc z (Ax^ = 70; 8a for 1 dof), allowing us to reject both hypotheses with high 
confidence. Interestingly, the "replenishing dust" model is nearly indistinguishable from 
an model because the dimming is directly proportional to distance traveled and thus 
mathematically quite similar to the effects of a cosmological constant. Consequently, we 
cannot discriminate this model from an fi^-dominated model strictly from its behavior in 
the magnitude-redshift plane (and probably never will be able to, given the small magnitude 
differences). However, the fine tuning required of this dust's opacity, replenishing rate, and 
velocity (> 1000 km s~^ for it to fill space uniformly without adding detectable dispersion) 
makes it unattractive as a simpler alternative to a cosmological constant. 

4.3. Exploring Dark Energy 

Despite the results of the last section which favor the dark-energy interpretation of 
SNe la, we avoid using this conclusion as a starting point for exploring the nature of dark 
energy. To do so would be to engage in "circular reasoning" or to incur more free parameters 
than our limited dataset can usefully constrain. Instead, we embark on a parallel study 
from the previous section. Here we use distance-independent information to justify the 
cosmological interpretation of SNe la and combine with other experiments to study dark 
energy. 

The potential for luminosity evolution of corrected SN la distances has been studied 
using a wide range of local host environments. No dependence of the distance measures 
on the host morphology, mean stellar age, radial distance from the center, dust content, 
or mean metallicity has been seen (Riess et al. 1998; Perlmutter et al. 1999; Hamuy et al 
2000). No differences in the inferred cosmology were seen by Sullivan et al. (2003) for SNc la 
in early-type hosts or late-type hosts at high redshifts. These studies limit morphology 
dependence of SN la distances to the 5% level. Detailed studies of distance-independent 
observables of SNe la such as their spectral energy distribution and temporal progression 
have also been employed as probes of evolution; see Riess (2000), Leibundgut (2001), 
and Perlmutter & Schmidt (2003) for reviews. The consensus interpretation is that there 
is no evidence for evolution with limits at or below the statistical constraints on the 
average high-redshift apparent brightness of SNe la. The observed nominal dispersion 
of high-redshift SNe la substantially limits the patchiness of uncorrected extinction, and 
near-IR observations of a high-redshift SN la demonstrate that a large opacity from grayish 
dust is unlikely (Riess et al. 2001). Non-SN constraints on gray dust from QSOs observed 
in X-rays (Paerels et al. 2002; Ninomiya, Yaqoob, & Khan 2003) and a partial or complete 



-28- 



resolution of the far-IR background by SCUBA (Chapman et al. 2003) place stringent 
hmits of less than a few percent of dimming at z — 0.5 from gray dust. 

Based on this evidence, we will adopt in the following analysis an a 'priori constraint 
that the net astrophysical contamination of SN la distance measures does not exceed their 
statistical uncertainty in their mean brightness. Quantitatively, our adopted limit on 
systematics is defined to be 5% per Az at z > 0.1. 

First we consider the SN data within an FRW cosmology of unknown curvature 
and mass density (with a flat prior on all parameters), with the simplest description 
of a dark-energy component (i.e., a cosmological constant) using equation (11). Joint 
confidence intervals in the Qm ~ plane were derived after numerical integration of 
the probability density P{Hq) oc cxp(— x;^(i/o)/2) over all values of Hq and are shown in 
Figure 8. Compared to the same analysis in Ricss ct al. (1998), the gold sample presented 
here reduces the area of the la contour by a factor of 6 (a factor of 7 including the silver 
sample). With the current sample, the 4cr confidence intervals (i.e., > 99.99% confidence) 
are now fully contained within the region where Qa > 0. The "concordance" model, of 
Qm = 0.27, = 0.73 lies within the la contour (though just outside of it with the addition 
of the silver set). For a flat geometry prior, we measure Q,m — 0.29±o;oi (equivalently 
Qa = 0.71). The iJST-discovered SNe la alone decrease the area of the la contour by a 
factor of 1.5 (in the gold sample) due to their high mean redshift. 

An alternative approach with good precedent (Garnavich et al. 1998; Perlmutter et al. 
1999) is to consider a flat Universe and a generalized dark-energy component parameterized 
by its (assumed) constant equation of state, w = P/ pc?. Flatness is assumed either on 
theoretical grounds (i.e., as a consequence of inflation) or on observational grounds from 
the characteristic angular size scale of the CMB fluctuations (Spergel et al. 2003, and 
references therein). In this case the luminosity distance is given by 

dL = CH^\1 + Z) r dz[{l + zf{QM) + (1 - ^m){1 + ^)3(l+-)]-V2. (13) 

Jo 

We determined the probability density in the flu — w plane in the same manner 
as above and the results are shown in the left panel of Figure 9. The SN la data alone 
require w < —0.5 for any value of Q^m at the 95% confldence level and are consistent with 
w — —1 (i.e., dark energy resembhng a cosmological constant) at the 68% confldence level 
for 0.20 < Qm < 0.35. 

Utilizing SN-independent constraints in this plane (primarily to constrain Qm) yields 
far more precise constraints for w due to the strong degeneracy between w and for 
SNe la. In the left panel of Figure 9 we use Qm = 0.27 ± 0.04 (at z = 0) as a simple 
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approximation to the constraints derived from numerous SN-independent experiments (see 
Frccdman & Turner 2003 for a review). Alternatively (right panel, Figure 9), we used the 
WMAPext (Bennett et al. 2003; Spergel et al. 2003) measurement of the reduced distance to 
the surface of last scattering ai z = 1089 and the Two-Degree Field Galaxy Redshift Survey 
(2dFGRS) measurement of the growth parameter, / = {Qm/{^m + (1 - f^Af)(l + -^)^"'))°•^ 
to derive independent constraints in the Qm — w plane. The results from either approach are 
similar. We find w — — 1.02±[] and w — — 1.08±^;f^ with the use of an Um prior and from 
WMAPext -|-2dFGRS, respectively. Using the prior on Qm, the 95% confidence interval is 
-0.78 >w> -1.46, or w < -0.76 ii w > -1. The 99% level constrains w < -0.72. These 
results are somewhat more constraining than those from Tonry et al. (2003), Barris et al. 
(2004), Knop et al. (2003). 

The high relative redshifts of the iJST-discovered SNe la provide little additional 
power to constrain a w parameter which is fixed a priori to be redshift-independent. The 

precision of such a constrained study of dark energy is most sensitive to the sheer number 
of SNe la with a relatively weak dependence on their redshifts (a useful approximation is 
Am ^ [1.82;/(1 + e^)]Aw for small Aw, w ~ —1, and z < 2) until the systematic error limit 
is reached. 

A more ambitious and potentially more revealing approach to studying dark energy 
is to allow for both an unconstrained value of the equation of state (at some fiducial 
redshift, e.g., z = 0) and its time evolution, i.e., w{z) = Wq +w'z, where w' = ^|^=q. 
This parameterization provides the minimum possible resolving power to distinguish a 
cosmological constant and a rolling scalar field from their time variation (or lack thereof). 
Indeed, rejection of the hypothesis that w' — would rule out a cosmological constant as 
the dark energy (as would the determination that w ^ —1). The measured value of w' 
would provide an estimate of the scale length of a dark-energy potential. The only previous 
estimate of w', by Di Pietro & Claeskens (2003), used the set of SNe la from Perlmutter et 
al. (1999) and the constraints Qtotai = 1 and Qm = 0.3, and concluded —12 < w' < 12 at 
the 95% confidence level (best fit: wq — — 1.4,w' = 2.3). 

For w{z) — Wo + w'z, we employ (following Linder 2003) 

dL = CH^\1 + Z) r dz[{l + zf{nM) + (1 - ^^m)(1 + ^)3(l+-o-')e3-'^]-l/2. (14) 

Jo 

A strong degeneracy exists among the three free parameters wq, w', and Qm, requiring the 
use of independent experimental constraints in this space to make progress. To this end 
we use the previous prior, Qm — 0.27 ± 0.04 (at z — 0), which has the advantage of being 
independent of redshift while providing a good approximation to all non-SN cosmological 
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constraints in this space. 

We have avoided using the CMB measurements directly as an additional constraint on 
the time evolution of the equation of state of dark energy due to difficulties which arise in the 
analysis of the CMB at 2; >> 1 where the linear expansion of the epoch-dependent equation 
of state diverges. Because constraints on w{z) from the CMB are derived from an integration 
between z — and z = 1089, diverging formulations of w{z) are unsuitable. Linder (2003) 
has proposed a more stable parameterization of w{z) = wq + Wa{z/{1 + z)); however, for 
large values of Wo + Wa, which are not rejected by the Di Pietro & Claeskens (2003) analysis, 
the CMB integral remains ill-behaved. In addition, the CMB measurements provide little 
direct leverage on Wa or w' compared to SNe la at 2; ~ 1 (because ^oarkEnergy ~ at 
z » 1). Therefore, we used the w' parameterization which is well-suited to our SN sample, 
and a simple prior on which avoids problems evaluating functions involving w{z) at 
z » 1. 

In Figure 10 we show constraints in the Wq — w' plane (after marginalizing over 0^)- 
Using the gold subset, we find Wq — — 1.31±q;|| and w' — 1.48±Q;gQ with the uncertainties 
in both parameters strongly correlated. A cosmological constant (i.e., wq — —1, w' — 0) is 
separated from the best fit along the direction of the major axis of the error ellipse, lying at 
the boundary of the joint 68% confidence level. If we constrain the recent behavior of dark 
energy to be like a cosmological constant (i.e., Wq = —1) we find w' = 0.60 ± 0.47. Models 
with Wo < —1 or "phantom energy" (Caldwell et al 2003) violate the Dominant-Energy 
condition {p + p > 0) and are extremely speculative at this point (and may rip apart the 
Universe in the future), so if we constrain the analysis to wq > —1 we find w' = 0.6 ± 0.5 
and Wq < —0.72 with 95% confidence. Unfortunately, few theoretical predictions about the 
size of w' exist for dynamic models of dark energy. However, we can reject the possibility 
that dark-energy evolution is currently very rapid (i.e., \w'\ > a few). This conclusion alone 
limits the rate at which simple rolling scalar fields could reach their true minima. The 
consequences of this statement for predicting the future fate of the Universe are discussed 
in §5. 

Our new constraints in the w' — w plane provide a substantial factor of ~ 8 improvement 
over the analysis by Di Pietro & Claeskens (2003) of the SCP (Perlmutter et al. 1999) 
data. Still, greater precision for the measurement of w' is needed before the proximity (or 
separation) from = would provide a compelling, empirical case for (or against) a static 
dark energy (i.e., a cosmological constant). The addition of the silver sample has only a 
modest impact on this analysis, as seen in Figure 10 (although a cosmological constant 
crosses to just outside the nominal 68% confidence interval). 

The iJ/^ST-discovered SNe la provide significant leverage in the Wq — w' plane due to 
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their high mean redshift. Figure 10 (upper left panel) shows the w' — w plane without 

inchiding the HST-discovered objects. The iKT-discovered SNe la alone increase the 
precision (i.e., reduce the area of the confidence intervals) by an impressive factor of 
1.9, although they account for only 10% of the sample. Previous studies in support of a 
dedicated, space-based mission to measure wq and w' have concluded that a SN la sample 
must extend to ^ > 1.5 to adequately break degeneracies in this parameter space (hinder 
& Huterer 2003), a conclusion supported by our analysis. The current relative dearth of 
SNe la at z > 1 compared to their number at z < 1 indicates that significant progress 
can still be made in the constraints on w' . Proposals for a Supernova Acceleration Probe 
(SNAP) or a Joint Dark Energy Mission (JDEM) predict an improved constraint for w' 
over our current analysis by a factor of 3 to 4, assuming a similar-sized improvement in our 
knowledge of VLm from the Planck Satellite and a continued ability to reduce systematic 
errors (hinder & Huterer 2003). However, the current sample is rapidly growing in size and 
we may expect progress in our constraints on the nature of dark energy in the next few 
years. 

5. Discussion 

5.1. Cosmological Constraints 

SNe la at z > 1 provide valuable and unique contributions to our current understanding 
of the cosmological model. The current sample of such SNe la, though greatly expanded 
here, remains small (i.e., < 10). Our capacity to constrain simultaneously the full debated 
range of cosmological and environmental parameters is therefore limited. Consequently, we 
have chosen to test specific and narrow questions in the context of well-defined assumptions 
or in conjunction with independent information. It is important to recognize that the 
conclusions garnered from any analysis cannot furnish a priori information for a subsequent 
analysis. Readers should carefully consider which priors they are using and where they 
came from before selecting which analysis presented here provides a relevant incremental 
gain. 

The two most extreme analyses (in the sense of the breadth of their priors) presented 
here also realize the most significant gains from the addition of our highest-redshift SNe la 
discovered with HST. The kinematic (i.e., cause- independent) interpretation of SN la 
distances and redshifts (independent of all other experiments) is most consistent with two 
distinct epochs of expansion: a recent accelerated expansion and a previous decelerated 
expansion with a transition between the two at z ~ 0.5. This is a generic requirement of 
a mixed dark matter and dark energy universe, and it may even be a feature of unrelated 
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cosmological paradigms (which are beyond our scope to consider here). The data are 
not consistent with many astrophysical interpretations posited in heu of dark energy. 
Notable examples include the attenuation produced in a universe filled with gray dust at 
2; > 1, or a luminosity evolution which is a simple, monotonic function of redshift. These 
interpretations are robust against the exclusion of any individual SN la used in the analyses 
and therefore represent an improvement over the results of Riess et al. (2001). 

A vacuum-driven metamorphosis model (VCDM) has been proposed by Parker & Raval 
(1999) to explain the cause of accelerated expansion. In this model the Universe makes a 
transition to a constant-scalar-curvaturc (induced by a quantized non-interacting scalar field 
of very small mass in its vacuum state) at 2; ~ 1. This model differs from a quintessence 
model in that the scalar field is free and thus interacts only with the gravitational field. 
The transition as described by Parker, Komp, & Vanzella (2003) is far more abrupt than for 
a cosmological constant and is characterized by Wq — —1.3 and w' — —0.8 for values of 
within its likely range. These values lay just outside the 99.5% confidence level. However, 
other future variations of the central idea of this model may yield better fits to the data. 

The "replenishing dust" model addressed in §4.2 is one example of a variety of 
astrophysical dimming models for which most of the apparent dimming occurs by 2; ~ 0.5, 
with little additional dimming at z > 0.5. Such models require an additional parameter or 
a logarithmic parameterization to dampen the dimming and fit the new SNe la at higher 
redshifts. Another model with this behavior would be evolution which is proportional to 
look-back time (Wright 2002). While possible, such dimming behavior, especially if in 
the form of luminosity evolution, would seem implausible. We may expect evolution (or 
dust production) to be coupled to the observed evolution of stellar populations, galaxy 
morphologies, sizes, large-scale structure, or even chemical enrichment. None of these 
known varieties of evolution are largely completed hy z — 0.5 starting from their properties 
at 2; = 0; quite the contrary, most of them have hardly begun, looking back to z — 0.5. 
A strong empirical argument against recent luminosity evolution is the independence of 
SN la distance measures on gross host morphology (Riess et al. 1998; Sullivan et al. 2003). 
The range of progenitor formation environments spanned by SNe la in early-type and 
late-type hosts greatly exceeds any evolution in the mean host properties between z = 
and z = 0.5. In the end, however, the only "proof" against astrophysical contamination 
of the cosmological signal from SNe la is to test the results against other experiments, 
independent of SNe la. 

SNe la are no longer unique in their requirement of a dominant dark-energy component. 
The bevy of SN-independent cosmological experiments, most notably the CMB, LSS and 
the ISW effect, provide a prior constraint of Qtotai — 1 and ~ 0.3. These priors preserve 
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the limited discriminating power of SN la data for resolving the nature of dark energy. 
What is its equation of state and has it been evolving? The constraints obtained here 
provide a substantial improvement in our ability to answer the latter (i.e., not rapidly), but 
the results are far from a compelling, empirical case for a cosmological constant or evolving 
dark energy. 

A parametric reconstruction of w{z) by Alam et al (2003) of the full set of SN la 
distances from Tonry et al (2003) and Barris et al (2003) concluded that dark energy 
"evolves rapidly" and has "metamorphosized" from w~0at2;~ltow<— lat2;~0. 
Their analysis differs from ours in a few important ways. Alam et al (2003) use a significant 
number of SNe that would fail our selection criteria of our gold sample including 26 SNe 
from the SCP (now excluded or flagged by Knop et al 2003) and another 13 from Barris et 
al (2003) of dubious reliability. Alam et al (2003) also use a different parameterization for 
w{z) and slightly different external cosmological constraints. These differences do account 
for a small part of the difference in our conclusions. However, despite these differences our 
results are fairly similar to theirs without the inclusion of the new H ST- discovered SNe la 
presented here as seen in the upper left panel of Figure 10. Without the HST-discovered 
SNe, the 95% confidence region resides in the quadrant in which Wq < —1 and w{z) was less 
negative in the past. Yet, with the addition of the HST-discovered SNe la, the contours 
shrink in area by a factor of 2 and shift in the direction of a cosmological constant, a 
position in the Wq — w' plane which is now at the boundary of the ~ 1 cr contour. Such 
large changes in the size and position of the error contours with the addition of new data 
indicates the crucial need for even more data before firmer conclusions can be reached. 

Knowledge of the equation of state of dark energy and its time evolution has profound 
implications for determining the fate of the Universe. For example, a joint constraint of 
Wq > —1 and < could provide the signature of a future recoUapse (i.e., "Big Crunch") 
for a linear potential field which becomes negative in the future (Kallosh & Linde 2003). 

Interestingly, the data indicate that this quadrant of the Wq — w' plane is the least favored. 
Values (in this quadrant) farthest from the origin (wq = —1 and w' = 0) would forecast 
the earliest possible recollapse and are the least favored by our data. Taking the Kallosh & 
Linde (2003) toy model at face value, our constraint would rule out such a future recollapse 
of the Universe in less than ~ 30 or 15 Gyr at the 95% or 99% confidence levels, respectively. 
This is somewhat more reassuring than the previous minimal remaining time to recollapse 
of 11 Gyr allowed by the observation that the dark-energy density has had time to evolve 
to Qa < 0.72 (Kallosh et al. 2003). 

A qualitatively different and more speculative fate, a hierarchical "ripping" of 
progressively smaller bound systems (i.e., "Big Rip") might occur if w{z) evolves to or 
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remains at values less than —1 (Caldwell, Kamionkowski, & Weinberg 2003). In this 
case the dark-energy density within any bound system would increase without limit as 
p oc a~^(^+^\ overcoming its binding energy. This eventuality may be foretold from joint 
constraints oi wq < —1 and > and (an assumed) linear extrapolation of the dark-energy 
scalar field. The current data are not inconsistent with this quadrant of the wq — w' plane 
(and will remain consistent if dark energy is a cosmological constant and measurements 
have finite precision). For an assumed constant equation of state, our 95% confidence 
interval would limit a "Big Rip" to be no sooner than 25 Gyr from now, only slightly more 
reassuring than the previous estimate of minimal remaining time to the Rip of 22 Gyr 
estimated by Caldwell, Kamionkowski, & Weinberg (2003). However, if w{z) is evolving to 
progressively more negative values (a result empirically consistent with our data), a Rip 
could occur far sooner. We are limited by our understanding of how the future equation of 
state of dark energy would track the scale factor, time, or some other global parameter of 
the Universe. Therefore, if dark energy is evolving in this direction, we cannot place any 
meaningful empirical limit on the minimum time to the Rip. 

We caution against naively extrapolating our empirical dynamical constraint (fit over 
a very small range in redshift) far from the time at which it was derived and is likely to be 
meaningful. The empirical approach to determining the fate of the Universe must rely on 
large extrapolations until dark energy is better understood. 

5.2. Tests of the Utility of Supernovae 

Many more SNe la at 2; > 1 may yield a more precise probe of dark energy, provided 
the statistical power of a larger sample can be realized. It is therefore of interest to use the 
current set of SNe la at 2; > 1 to consider whether there are any reasons to believe this goal 
could not be accomplished. 

The dispersion of the SNe la at 2; > 1.0 around the best fit model is 0.29 mag, similar 
to the 0.27 mag dispersion of the sample at 0.1 < z < 1.0. (We note that these values are 
higher than the nominal intrinsic dispersion of ~0.15 mag due to the sparse sampling and 
noisy photometry of SNe la observed at high redshifts.) In addition, the average per SN 
at 2; > 1.0 (9 objects in the gold sample) is 0.95 for the concordance model. We conclude 
that there is no evidence of excess dispersion for SNe la at 2; > 1. 

Another test of evolutionary effects is to compare the distribution of light-curve shapes 
for the //ST-discovered SNe la with those at low redshift. In Figure 11 we compare the 
individually normalized (i.e., after subtraction of the individual peak day and apparent 
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magnitude) , rest- frame (i.e., after K-corrections and correcting for 1 + z time dilation) 
light curves of these samples. Statistically, the measured magnitudes of iiZ/S'T-discovered SN 
light curves (with median redshift 1) are consistent with having been drawn from a parent 
population at low redshift. The fitted distributions of light-curve shape parameters from 
the MLCS2k2 parameterization, A (Riess et al. 1996a; Jha et al. 2004a), are also consistent 
for the two populations as shown in Figure 12. 

Some evidence for anomalously blue colors of the high-redshift SNe la in the objects 
from Riess ct al. (1998) was noted by Falco et al. (1999) and Leibundgut (2001). The 
indicated mean difference was a few hundredths of a magnitude in data transformed to 
rest-frame B — V with a significance approaching 2a. No such difference was seen in the data 
from Perlmutter et al. (1999), although the individual precision of their measured colors 
may have precluded the detection of such an anomaly. Here we have used the expanded 
data set and its redshift leverage to look for evidence of color evolution in both rest-frame 
B — V and U — B. The ultraviolet colors of SNe la are predicted by theory to be the most 
sensitive to chemical abundances (Hoffich, Wheeler, & Thielemann 1998). Figure 13 shows 
both of these colors as a function of redshift for the sample in a modest range around the 
unreddened colors of SNe la. In Figure 14 we have constructed histograms for the colors in 
three redshift ranges: z < 0.1, 0.1 < 2; < 0.6, and 0.6 < 2; < 1.6. For B — V, the middle 
bin is bluer in the mean by 0.02 mag than the lowest bin, and despite the increased sample 
size the significance of this difference remains 2a. This result is not surprising because 
it makes use of similar data as the previous analysis. Interestingly, the size of the mean 
differences decreases for the highest-redshift bin to 0.01 mag with an even lower significance 
of la. For the U — B data, the lowest and highest-redshift bins are consistent in their 
means of —0.45 ± 0.02 mag and —0.47 ± 0.02 mag, respectively. (Negligible data ioi U — B 
in the middle redshift bin is available.) Of the three comparisons presented here, only the 
previously noted B — V comparison between low and intermediate redshifts is suspect. If 
the colors of SNe la were evolving, we would expect the comparison at higher redshifts, and 
additionally with more chemically sensitive colors, to provide amplified evidence for the 
effect. In addition, the increased sample size in the low and intermediate-redshift bins has 
not produced an increase in the significance of the anomaly as would be expected for a real 
effect. We therefore conclude there is insufficient evidence for color evolution of SNe la with 
the current sample, but we encourage future investigators to remain aware of (and test for) 
this possibility. 

Overall, the data presented here suggest that SNe la at 2; > 1 remain empirically 
well-behaved and show promise for providing robust cosmological measurements. 

To test the sensitivity of our results to the light-curve fitting method used, we 
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compared the MLCS2k2 distances of the iKT-discovered SNe la to those derived from 
BATM (Tonry et al. 2003). For this comparison, the distance scales of the two methods 
were normahzed by matching the mean distances calculated for the same SNe la from 
Tonry et al (2003). In Figure 15 we show a comparison for the SNe la fit by both methods 
over a wide range in redshift. The mean difference (weighted by the quadrature sum of the 
individual uncertainties) for the if^T-discovered SNe la is 0.047 mag (in the sense that the 
BATM distances are larger). This difference is consistent (at the la confidence level) with 
an independent mean uncertainty of 0.18 mag for each method. If we limit the comparison 
to the 9 iJST-discovered SNe la at z > 0.9, the mean difference diminishes to < 0.01 
mag. Because this higher-redshift range provides the leverage utihzed for the updated 
cosmological conclusions presented here, we conclude that our results are insensitive to the 
light-curve fitting method. 

The problem of heterogeneity of SN la data may be addressed by future, massive SN 
surveys such as ESSENCE (Smith et al. 2002; Garnavich et al. 2002) or the CFHT Legacy 
Survey (Pain et al. 2002), which will obtain large sets of SNe la over the full range of 
explored redshifts (obviating the need for older, heterogeneous data). Until then, we must 
rely on a careful and judicious compilation of the available data and a clear description of 
how they were compiled. 

Farrah et al. have used the statistical measurement (most individual hosts were not 
detected) of sub-millimetre emission from high-redshift SN la hosts to infer the average 
internal extinction of these galaxies. They then assume that the line-of-sight extinction 
of supernovae matches the inferred mean extinction of the hosts. They conclude that the 
SNe la from hosts at higher redshifts may suffer more extinction than their lower redshift 
counterparts and this effect may bias the inferred cosmology. However, it has not been 
demonstrated that radio emission correlates well (if at all) with SN extinction. In contrast, 
apparent reddening has been shown to provide a significant reduction in distance dispersion 
and should account for any change in the relative proportions of host extinction along the 
line-of-sight. More relevant to this work, a trend of increased and uncorrected extinction 
with redshift does not match the apparent magnitudes of the new SNe la presented here in 
the highest redshift bin. 

5.3. Lensing 

Until now, we have evaluated the SN data in the context of a homogeneous model 
for the matter distribution, also known as the "filled-beam" model. The inhomogeneity of 
matter along the SN line-of-sight distorts space-time, directing more or less photon-filled 
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lines-of-sight to our detectors depending on the matter content along the hne of sight. 
Averaged over enough hnes-of-sight. SN flux is conserved, and the inferred cosmology is 
unaffected. However, further investigation of possible lensing is warranted for small sample 
sizes and in regard to possible selection biases. 

We have estimated the expected lensing along the line-of-sight of each of the 
ifST-discovered SNe la as well as for 100 randomly selected positions in each of the GOODS 
fields. The lensing is estimated using the same multiple lens-plane methodology employed 
by Benftez et al. (2002) for SN 1997ff. Photometric and spectroscopic rcdshifts of all 
foreground lenses were derived from the GOODS catalogs. Their masses were estimated 
from the rest-frame 5-band luminosities and from the TuUy-Fisher and Faber- Jackson 
relations corrected for evolution (Ziegler et al. 2001; Boehm et al. 2003; Kochanek 1996; 
Treu et al. 2002). Because SN luminosity distances here and elsewhere have been modeled 
using "fiUed-beam" cosmological models, we have calculated the amphfications relative to a 
fiUed-beam model and for = 0.3, JIa = 0.7. The results are shown in Figure 16. 

As expected, the distribution of net amplifications for random 3-dimensional positions 
in the fields scatters about unity with a width steadily increasing with redshift (simulations 
predict a — 0.03 mag per unit redshift as indicated; D. Holz 2003, private communication). 
Relatively underdense and frequent lines-of-sight resulting in net deamplification are 
compensated by an extended tail of strong magnifications, resulting in an increasingly 
skewed distribution at higher redshifts. These results are in good accord with Monte Garlo 
simulations of this effect (Holz 1998; Metcalf & Silk 1999). 

We find that the magnification estimates for the SNe la discovered in the GOODS 
survey are consistent with having been drawn from those in the random positions sample. 
We conclude that there is no evidence of a lensing-induced selection bias of our sample. 
This is not unexpected since most of our SNe la were discovered at > 2 mag brighter than 
the survey limit and lensing would contribute an insignificant < 0.2 mag of amplification 
for 98% of the population at z < 1.5. 

Interestingly, there is an indication that the discovery of SNe la in the HDF-N found 
by other surveys could have been favored by amplification. The net amplification of SN 
1997ff of ~ 0.3 mag (Benitze et al. 2002) is unusually large and the SN was within 0.5 mag 
of its survey hmit (Gilliland, Nugent, & Phillips 1999). It is less clear whether the discovery 
of SN 2002dd {z — 0.95) could have been favored by its 0.2 mag amphfication due to the 
serendipity of its discovery (Blakeslee et al. 2003). 

In the absence of an apparent selection bias (or bad luck) in our SN survey, it 
is unnecessary (and even undesirable) to correct individual SNe for the estimated net 
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amplification, as the sample will approach a mean amplification of unity when considered 
in a filled-beam model. An important exception would be to correct (or flag) SNe whose 
predicted amplification places them on the strong lensing tail, such as SN 1997ff as done 
by Benitez et al (2001), Tonry et al (2003) and here. However, an independent test of this 
conclusion is available by comparing the predicted amplifications with the observed residuals 
from a good-fit cosmological model. In Figure 17 we show this comparison and their 
derived correlation. Empirically we find the residuals to be 1.6 ± 0.9 times the predicted 
magnifications, consistent with the theoretical relation (i.e., unity) and inconsistent with 
no correlation at the 1.7a (80%) confidence level. Much of the leverage comes from SN 
1997ff; without it, the correlation is 1.0 ± 1.2. As the sample expands, it may soon be 
possible to provide independent evidence of the correlation of dark matter and light from 
the high-redshift SN sample. 

Mortsell, Gunnarsson, & Goobar (2001) (hereafter MGG2001) and more recently 
Gunnarsson (2004) have made predictions of the expected magnification of two of the SNe la 
in our sample, both contained in the original HDF-N. For SN 1997ff {z = 1.7) these authors 
concluded the magnification was uncertain and potentially large (i.e., as high as a factor 
of a few). In contrast, Benitez et al. (2002) estimated a much more modest magnification 
of 0.34 ± 0.12 mag, in good agreement with the values estimated by Lewis & Ibata (2001) 
and Riess et al. (2001). The source of the potentially higher magnification can be traced 
to differences in the treatment of the mass scahng of the foreground lenses. MGG2001 
treated the lensing galaxies as unevolved with an unknown but possibly high mass scaling 
(e.g., velocity dispersions of 200 to 300 km s~^ for galaxies with Mb — —19.5 + 51og/i). 
Benitez et al. (2002) used the B-band Tully-Fisher relation observed by Ziegler et al. 
(2002) for late-type galaxies at redshifts 0.1-1.0 to provide a slope and normalization which 
accounts for evolution (which yields a velocity dispersion M^, ~ 150 km s^^). The result is 
a much smaller mass scale than that considered by MGG2001 for the foreground lenses of 
high-redshift SNe. 

For SN 2003es, Gunnarsson (2004) estimated a magnification factor of 1.15 (with a 
velocity dispersion of 170 km s^^) for 13 apparent foreground galaxies (at z < 0.968) based 
primarily on optical photometric redshifts. To span the rest-frame optical breaks in the 
spectral energy distributions of potential high-redshift foreground galaxies and reliably 
estimate their redshifts, it is important to use near-IR observations with good precision. 
We reanalyzed the photometric redshifts of these 13 galaxies using the high-resolution 
near-IR data from NICMOS (Budavari et al. 2000). In the majority of cases we found that 
the data from Fernandez-Soto, Lanzetta & Yahil (1999) and Gwynn & Hartwick (1996) 
contained detections in only 3 or 4 bands whereas the NICMOS data provided a total of 5 
to 7 bands and with greater leverage. We found 7 of the 13 galaxies to lie in the background 
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of the SN. Repeating the analysis of Gunnarsson (2004) using their Q-LET algorithm for 
the 7 remaining foreground galaxies reduced the predicted magnification factor to 1.10. 
However, calculations from the Q-LET algorithm are made relative to a non-filled-beam 
cosmology and hence all lines-of-sight are amplified (in the absence of any foreground lenses 
the (5-Li?r amplification would be 1.0). At the redshift of SN 2003cs, the average strong 
lensing amplification would be ~0.03 mag so the remaining excess would be a factor of 
1.07. The remaining difference between this value and our estimate of 1.03 results from the 
previously discussed difference in the treatment of the size and evolution of the mass scale. 
Our analysis, comparing the predicted magnification and observed cosmological residuals, 
is consistent with the mass scahng from Ziegler et al. (2002) as utihzed here (and would be 
inconsistent with a scale approximately twice as large). 

It is tempting to consider that we have reached the end of the beginning in the 
exploration of dark energy. Two reliable and independent routes require it in addition to 
a third more tentative investigation via the integrated Sachs- Wolfe effect (Scranton et al. 
2003). SNc la continue to provide the most direct route to illuminating dark energy because 
their light can be measured propagating from within its era of dominance. Two clues about 
dark energy, its equation of state and its recent time evolution, would be invaluable aids to 
test and provoke theories. We suggest that the most efficient way forward in the near term 
is by simultaneously mining both ends of the observable redshift range: at 2; < 1 generally 
from the ground, and at 2; > 1 generally from space. The constraints presented here in the 
Wq — w' plane have reduced the allowable range of w' from a factor of ~ 10 to ^ 1 while 
retaining the constraints on wq within —1.4 < wq < —0.7. With continued determination, 
an improvement in precision by a factor of a few in this plane is expected. 



6. Summary and Conclusions 

We have conducted the first space-based SN search and follow-up campaign using the 
ACS on board HST. The search parameters and the full list of 42 new SNe are provided 
elsewhere (Strolger et al. 2004). We reviewed the sample of SNe la harvested from the 
survey and examined its cosmological significance. The key results can be summarized as 
follows. 

(1) We obtained multi-color light curves and spectroscopic redshifts for 16 new SNe la 
which uniformly sample the redshift range 0.2 < z < 1.6. Twelve of these are classified by 
their spectra, 2 from their red, early-type host galaxies, and 2 by photometric diagnostics. 
Three of the SN spectra are at the highest redshifts yet observed for SNe. Six of the SNe la 
are among the seven highest- redshift known; all are at z > 1.25. These data provide a 
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robust extension of the Hubble diagram to 1 < z < 1.6. 

(2) Utilizing a simple kinematic description of the magnitude-redshift data, we find 
that the SNe la favor recent acceleration and past deceleration at the 99.2% confidence 
level. An alternate kinematic parameterization requires a positive jerk (third derivative of 
the scale factor) . The best- fit redshift of the transition between these kinematic phases is 
z — 0.46 ± 0.13, although the precise value depends on the kinematic model employed. 

(3) We have compared the goodness-of-fit of cosmological models and simple models of 
astrophysical dimming. The "gold" sample of 157 SNe la is consistent with the "cosmic 
concordance" model (Qjv^ = 0.3, Qa = 0.7) with x^^jj = 1.06. The data reject at high- 
confidence simple, monotonic models of astrophysical dimming which are tuned to mimic 
the evidence for acceleration at z ~ 0.5. These models include either a universe filled with 
gray dust at high redshift, or luminosity evolution oc z. More complex parameterizations of 
astrophysical dimming which peak at 2; ~ 0.5 and dissipate at 2; > 1 remain consistent with 
the SN data (but appear unattractive on other grounds). 

(4) For a flat Universe with a cosmological constant, we measure = 0.29±q;q3 
(equivalently, = 0.71). When combined with external flat-Universe constraints including 
the CMB and LSS, we flnd for the dark-energy equation of state w = —1.02±q \1 (and 

w < —0.76 at the 95% confidence level) for an assumed static equation of state of dark 
energy, P — wpc^. 

(5) Joint constraints on both the recent equation of state of dark energy and its time 
evolution are a factor of ~ 8 more precise than its first estimate and twice more precise 
than those derived without the SNe la discovered by HST. Both of these dark energy 
properties are consistent with a cosmological constant (i.e., with wq — —1.0, w' — 0) and 
are inconsistent with very rapid evolution of dark energy (i.e., \w'\ > a few). The absence 
of rapid evolution places constraints on the time in which a simple scalar field could evolve 
to recollapse the Universe. Specifically, the timescale to a potential recoUapse is larger 
than ~30 Gyr. If dark energy is evolving towards more negative w, we cannot place any 
meaningful limit on the minimum time to a (speculative) Big Rip. 
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A. Appendix: The Full Sample 

Distance measurements to individual SNe depend on the algorithms used to estimate 
their K-corrections, fit their light curves, and infer their extinction. There is currently 
no single set of algorithms which are considered by consensus to be the optimal ones. 
Rather, different methods may have advantages depending on the breadth and quality of 
the observational record of any individual SN la. In addition, algorithms improve as their 
training samples grow. Here we present the full cosmological sample of SNe la used in this 
work in Table 5. Their virtue is that all distance estimates were derived from a single set of 
algorithms, MLCS2k2 (Jha et al. 2004a), with the broadest set of training data available at 
this time including f/-band data. Additional advantages include a consistent and thorough 
reanalysis of quality criteria for all currently published SNe la, an exercise resulting in the 
rejection of many SNe from our "gold" sample whose observational records have one or 
more shortcomings (see §3.1 for discussion). 

The zeropoint, distance scale, absolute magnitude of the fiducial SN la or Hubble 
constant derived from Table 5 are all closely related (or even equivalent) quantities which 
were arbitrarily set for the sample presented here. Their correct value is not relevant for 
the analyses presented which only make use of differences between SN la magnitudes. Thus 
the analysis are independent of the aforementioned normalization parameters. 
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Table 4. Comparison of Gold Set Data to Models 



Model 


X2(for 157 SNe la) 


Qm = 0.27, Qa = 0.73 


178 


= 1.00, = 0.00 


325 


Qj^^ = 0.00, Qa = 0.00 


192 


High-redshift gray dust (with Qm = 1-00, Qa = 0.00) 


307 


Replenishing dust (with Qm = 1-00, Qa = 0.00) 


175 


Dimming oc z (with O^f = 1-00, Qa = 0.00) 


253 



Table 5. MLCS2k2 Full Sample 



SN 


z 


MS 


a* 


host Ay 


sample 





SN 1990T 


0.0400 


36.38 


0.19 


0.37 


gold 


SN 1990af 


0.050 


36.84 


0.21 


-0.04 


gold 


SN 1990O 


0.0307 


35.90 


0.20 


0.11 


gold 


SN 1991S 


0.0560 


37.31 


0.18 


0.20 


gold 


SN 1991U 


0.0331 


35.54 


0.20 


0.37 


gold 


SN 1991ag 


0.0141 


34.13 


0.25 


0.12 


gold 


SN 1992J 


0.0460 


36.35 


0.21 


0.25 


gold 


SN 1992P 


0.0265 


35.64 


0.20 


0.17 


gold 


SN 1992aq 


0.101 


38.73 


0.20 


-0.03 


gold 


SN 1992ae 


0.075 


37.77 


0.19 


0.16 


gold 


SN 1992au 


0.061 


37.30 


0.22 


0.09 


gold 


SN 1992al 


0.0141 


34.12 


0.25 


0.05 


gold 


SN 1992ag 


0.0262 


35.06 


0.24 


0.54 


gold 


SN 1992bl 


0.0430 


36.53 


0.19 


-0.04 


gold 


SN 1992bh 


0.0450 


36.97 


0.18 


0.35 


gold 


SN 1992bg 


0.036 


36.17 


0.19 


0.21 


gold 


SN 1992bk 


0.058 


37.13 


0.19 


0.03 


gold 


SN 1992bs 


0.063 


37.67 


0.19 


0.26 


gold 


SN 1992bc 


0.0186 


34.96 


0.22 


-0.04 


gold 


SN 1992bp 


0.079 


37.94 


0.18 


0.03 


gold 


SN 1992br 


0.088 


38.07 


0.28 


-0.04 


gold 


SN 1992bo 


0.0178 


34.70 


0.23 


-0.01 


gold 


SN 1993B 


0.071 


37.78 


0.19 


0.36 


gold 


SN 1993H 


0.0251 


35.09 


0.21 


0.05 


gold 


SN 19930 


0.052 


37.16 


0.18 


0.13 


gold 


SN 1993ah 


0.0286 


35.53 


0.21 


0.26 


gold 


SN 1993ac 


0.0490 


36.90 


0.20 


0.54 


gold 


SN 1993ag 


0.050 


37.08 


0.19 


0.28 


gold 


SN 1993ae 


0.0180 


34.29 


0.23 


0.00 


gold 


SN 1994B 


0.089 


38.50 


0.17 


0.00 


silver 


SN 1994C 


0.051 


36.67 


0.16 


0.00 


silver 


SN 1994M 


0.0244 


35.09 


0.20 


0.23 


gold 


SN 1994Q 


0.0290 


35.70 


0.19 


0.33 


gold 


SN 1994S 


0.0161 


34.50 


0.24 


0.06 


gold 


SN 1994T 


0.0360 


36.01 


0.20 


0.09 


gold 


SN 1995E 


0.0116 


32.96 


0.29 


2.48 


silver 


SN 1995K 


0.478 


42.48 


0.23 


0.04 


gold 


SN 1995M 


0.053 


37.17 


0.15 


0.00 


silver 


SN 1995ap 


0.230 


40.44 


0.46 


0.00 


silver 


SN 1995ao 


0.300 


40.76 


0.60 


0.00 


silver 


SN 1995ae 


0.067 


37.54 


0.34 


0.00 


silver 
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SN 


z 




a* 


host Ay 


sample 


SN 1995az 


0.450 


42.13 


0.21 


— 


gold 


SN 1995ay 


0.480 


42.37 


0.20 


— 


gold 


SN 1995ax 


0.615 


42.85 


0.23 


— 


gold 


SN 1995aw 


0.400 


42.04 


0.19 


— 


gold 


SN 1995as 


0.498 


43.21 


0.24 


— 


silver 


SN 1995ar 


0.465 


42.81 


0.22 


— 


silver 


SN 1995ac 


0.0490 


36.52 


0.20 


0.40 


gold 


SN 1995ak 


0.0219 


34.70 


0.22 


0.56 


gold 


SN 1995ba 


0.3880 


42.07 


0.19 


— 


gold 


SN 1995bd 


0.0152 


34.11 


0.25 


0.70 


gold 


SN 1996C 


0.0276 


35.90 


0.20 


0.34 


gold 


SN 1996E 


0.425 


41.70 


0.40 


0.35 


gold 


SN 1996H 


0.620 


43.11 


0.30 


0.09 


gold 


SN 19961 


0.570 


42.81 


0.25 


0.14 


gold 


SN 1996J 


0.300 


41.01 


0.25 


0.23 


gold 


SN 1996K 


0.380 


42.02 


0.22 


0.02 


gold 


SN 1996R 


0.160 


39.08 


0.40 


0.00 


silver 


SN 1996T 


0.240 


40.68 


0.43 


0.00 


silver 


SN 1996U 


0.430 


42.33 


0.34 


0.08 


gold 


SN 1996V 


0.0247 


35.33 


0.25 


0.00 


silver 


SN 1996ab 


0.124 


39.20 


0.22 


0.00 


gold 


SN 1996bo 


0.0165 


33.82 


0.27 


0.77 


gold 


SN 1996bv 


0.0167 


34.21 


0.23 


0.71 


gold 


SN 1996bl 


0.0348 


36.17 


0.19 


0.33 


gold 


SN 1996cg 


0.490 


42.58 


0.19 


0.63 


silver 


SN 1996cm 


0.450 


42.58 


0.19 


— 


silver 


SN 1996cl 


0.828 


43.96 


0.46 


— 


gold 


SN 1996ci 


0.495 


42.25 


0.19 


— 


gold 


SN 1996cf 


0.570 


42.77 


0.19 


— 


silver 


SN 1997E 


0.0132 


34.02 


0.26 


0.12 


gold 


SN 1997F 


0.580 


43.04 


0.21 


— 


gold 


SN 1997H 


0.526 


42.56 


0.18 


0.45 


gold 


SN 19971 


0.172 


39.79 


0.18 




gold 


SN 1997N 


0.180 


39.98 


0.18 




gold 


SN 1997P 


0.472 


42.46 


0.19 




gold 


SN 1997Q 


0.430 


41.99 


0.18 




gold 


SN 1997R 


0.657 


43.27 


0.20 




gold 


SN 1997Y 


0.0166 


34.54 


0.23 


0.25 


gold 


SN 1997ai 


0.450 


42.10 


0.23 




gold 


SN 1997ac 


0.320 


41.45 


0.18 




gold 


SN 1997aj 


0.581 


42.63 


0.19 




gold 
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SN 


z 




(7* 


host Ay 


sample 


SN 1997aw 


0.440 


42.57 


0.40 


0.80 


gold 


SN 1997as 


0.508 


41.64 


0.35 


0.85 


gold 


SN 1997am 


0.416 


42.10 


0.19 


0.00 


gold 


SN 1997ap 


0.830 


43.85 


0.19 


— 


gold 


SN 1997af 


0.579 


42.86 


0.19 


— 


gold 


SN 1997bh 


0.420 


41.76 


0.23 


0.60 


gold 


SN 1997bb 


0.518 


42.83 


0.30 


0.11 


gold 


SN 1997bj 


0.334 


40.92 


0.30 


0.34 


gold 


SN 1997ck 


0.970 


44.13 


0.38 


0.17 


silver 


SN 1997cn 


0.0175 


34.52 


0.25 


0.01 


gold 


SN 1997cj 


0.500 


42.74 


0.20 


0.15 


gold 


SN 1997ce 


0.440 


42.08 


0.19 


0.08 


gold 


SN 1997dg 


0.0297 


36.12 


0.20 


0.28 


gold 


SN 1997do 


0.0104 


33.73 


0.33 


0.44 


gold 


SN 1997ez 


0.778 


43.81 


0.35 


— 


gold 


SN 1997ek 


0.860 


44.03 


0.30 


— 


gold 


SN 1997eq 


0.538 


42.66 


0.18 


— 


gold 


SN 1997ff 


1.755 


45.53 


0.35 


0.00 


gold 


SN 19981 


0.886 


42.91 


0.81 


0.95 


gold 


SN 1998J 


0.828 


43.61 


0.61 


0.49 


gold 


SN 1998M 


0.630 


42.62 


0.24 


0.75 


gold 


SN 1998V 


0.0170 


34.47 


0.23 


0.28 


gold 


SN 1998ac 


0.460 


41.83 


0.40 


0.48 


gold 


SN 1998ay 


0.638 


43.30 


0.36 


— 


silver 


SN 1998bi 


0.740 


43.35 


0.30 


— 


gold 


SN 1998be 


0.644 


42.78 


0.26 


— 


silver 


SN 1998ba 


0.430 


42.36 


0.25 


— 


gold 


SN 1998bp 


0.0104 


33.21 


0.32 


0.19 


gold 


SN 1998CO 


0.0171 


34.68 


0.24 


0.20 


gold 


SN 1998CS 


0.0327 


36.08 


0.19 


-0.03 


gold 


SN 1998dx 


0.053 


36.97 


0.18 


0.04 


gold 


SN 1998ef 


0.0170 


34.18 


0.23 


0.07 


gold 


SN 1998cg 


0.0234 


35.36 


0.20 


0.29 


gold 


SN 1999Q 


0.460 


42.56 


0.27 


0.23 


gold 


SN 1999U 


0.500 


42.75 


0.19 


0.04 


gold 


SN 1999X 


0.0257 


35.41 


0.20 


0.31 


gold 


SN 1999aa 


0.0157 


34.58 


0.24 


0.02 


gold 


SN 1999CC 


0.0316 


35.85 


0.19 


0.09 


gold 


SN 1999cp 


0.0104 


33.56 


0.31 


0.08 


gold 


SN 1999da 


0.0121 


34.05 


0.32 


0.58 


silver 


SN 1999dk 


0.0141 


34.43 


0.26 


0.20 


gold 
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SN 


z 




cr* 


host Ay 


sample 


SN 1999dq 


0.0136 


33.73 


0.26 


0.43 


gold 


SN 1999ef 


0.0380 


36.67 


0.18 


0.05 


gold 


SN 1999fw 


0.278 


41.00 


0.41 


0.26 


gold 


SN 1999fk 


1.056 


44.25 


0.23 


0.19 


gold 


SN 1999fin 


0.949 


43.99 


0.25 


0.11 


gold 


SN 1999fj 


0.815 


43.76 


0.33 


0.23 


gold 


SN 1999ff 


0.455 


42.29 


0.28 


0.19 


gold 


SN 1999fv 


1.19 


44.19 


0.34 


0.24 


gold 


SN 1999fh 


0.369 


41.62 


0.31 


0.70 


silver 


SN 1999£n 


0.477 


42.38 


0.21 


0.15 


gold 


SN 1999gp 


0.0260 


35.62 


0.20 


0.18 


gold 


SN 2000B 


0.0193 


34.59 


0.23 


0.28 


gold 


SN 2000bk 


0.0266 


35.36 


0.21 


0.19 


gold 


SN 2000cf 


0.0360 


36.39 


0.18 


0.21 


gold 


SN 2000CI1 


0.0233 


35.14 


0.21 


0.08 


gold 


SN 2000ce 


0.0164 


34.47 


0.23 


1.02 


silver 


SN 2000dk 


0.0164 


34.41 


0.24 


-0.05 


gold 


SN 2000dz 


0.500 


42.75 


0.24 


0.09 


gold 


SN 2000eh 


0.490 


42.41 


0.25 


0.20 


gold 


SN 2000ee 


0.470 


42.74 


0.23 


0.13 


gold 


SN 2000eg 


0.540 


41.96 


0.41 


0.12 


gold 


SN 2000ea 


0.420 


40.79 


0.32 


1.05 


silver 


SN 2000ec 


0.470 


42.77 


0.21 


0.13 


gold 


SN 2000fr 


0.543 


42.68 


0.19 


— 


gold 


SN 2000fa 


0.0218 


35.06 


0.21 


0.44 


gold 


SN 2001V 


0.0162 


34.13 


0.23 


0.28 


gold 


SN 200 Ifs 


0.873 


43.75 


0.38 


0.64 


gold 


SN 2001fo 


0.771 


43.12 


0.17 


0.05 


gold 


SN 2001hy 


0.811 


43.97 


0.35 


0.03 


gold 


SN 2001hx 


0.798 


43.88 


0.31 


0.31 


gold 


SN 2001hs 


0.832 


43.55 


0.29 


0.10 


gold 


SN 2001hu 


0.882 


43.90 


0.30 


0.12 


gold 


SN 2001iw 


0.340 


40.71 


0.27 


0.73 


gold 


SN 2001iv 


0.397 


40.89 


0.30 


0.91 


gold 


SN 200 liy 


0.570 


42.88 


0.31 


-0.04 


gold 


SN 2001ix 


0.710 


43.05 


0.32 


0.53 


gold 


SN 2001jp 


0.528 


42.77 


0.25 


0.10 


gold 


SN 2001jli 


0.884 


44.23 


0.19 


-0.01 


gold 


SN 2001jb 


0.698 


43.33 


0.32 


0.15 


silver 


SN 2001jf 


0.815 


44.09 


0.28 


0.23 


gold 


SN 2001jni 


0.977 


43.91 


0.26 


0.18 


gold 
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Figure Captions 

Figure 1: Discovery- image sections from ACS FS50LP images around each SN. Panels 
on the left and middle show the discovery epoch and the preceding (template) epoch, 
respectively. The panels on the right show the results of the subtraction (discovery epoch 
minus template). Arrows indicate position of the SNe. Image scales and orientations are 
given. 

Figure 2: Multi-color light curves of SNe la. For each SN la, multi-color photometry 
transferred to rest-frame passbands is plotted. The individual, best-fit MLCS2k2 model is 
shown as a solid line, with a ±lcr model uncertainty, derived from the model covariance 
matrix, above and below the best fit. 

Figure 3: Identification spectra (in fA) of 12 of the new iJST-discovered high-redshift 
SNe la, shown in the rest frame. Classification features are analyzed in §2.3. The data are 
compared to nearby SN la spectra of the same age as determined by the light curves (see 
Table 3). Classification of the 5 SNe without spectra (SN 20031v, z = 0.94; SN 2002fx, 
z = 1.40; SN 2002hp, z = 1.31; SN 2003ak, z = 1.55; and SN 2003aj, z = 1.31) are discussed 
in section §2.3 and §3.1. 

Figure 4: MLCS2k2 SN la Hubble diagram. SNe la from ground-based discoveries in 
the gold sample are shown as diamonds, i/^T-discovered SNc la arc shown as filled symbols. 
Overplotted is the best fit for a fiat cosmology: — 0.29, fl\ = 0.71. 

Figure 5: Left panels: Joint confidence intervals for a two-parameter model of the 
expansion history, q{z) = qo + zdq/dz, from SNe la. The upper left shows the constraints 
derived from the gold sample, the lower left includes both gold and silver samples. For either 
set, the data favor the quadrant with recent acceleration (go < 0) and past deceleration 
{dq/dz > 0) with high confidence. Lines of fixed transition redshift {q{zt) = 0) are shown. 
Panels on the right illustrate the likelihood function for the transition redshift derived from 
the same samples. 

Figure 6: Kinematic SN la residual Hubble diagram. Upper panel: SNe la from 
ground-based discoveries in the gold sample are shown as diamonds, if^T-discovered SNe la 
are shown as filled symbols. Bottom panel: weighted averages in fixed redshift bins are 
given for illustrative purposes only. Data and kinematic models of the expansion history 
are shown relative to an eternally coasting model, q{z) — 0. Models representing specific 
kinematic scenarios (e.g., "constant acceleration") are illustrated. 

Figure 7: SN la residual Hubble diagram comparing cosmological models and models 
for astrophysical dimming. Upper panel: SNe la from ground-based discoveries in the gold 
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SN 


z 


MS 


a* 


host Ay 


sample 


SN 2001kd 


0.935 


43.99 


0.38 


0.14 


silver 


SN 2002P 


0.719 


43.22 


0.26 


0.11 


silver 


SN 2002ab 


0.422 


42.02 


0.17 


0.10 


silver 


SN 2002ad 


0.514 


42.39 


0.27 


0.09 


silver 


SN 2002dc 


0.475 


42.14 


0.19 


0.23 


gold 


SN 2002dd 


0.95 


44.06 


0.26 


0.24 


gold 


SN 2002fw 


1.30 


45.27 


0.19 


0.21 


gold 


SN 2002£x 


1.40 


45.09 


0.45 


0.49 


silver 


SN 2002hr 


0.526 


43.01 


0.27 


0.74 


gold 


SN 2002hp 


1.305 


44.70 


0.22 


0.19 


gold 


SN 2002kc 


0.216 


40.33 


0.18 


1.29 


silver 


SN 2002kd 


0.735 


43.09 


0.19 


0.21 


gold 


SN 2002ki 


1.140 


44.84 


0.30 


0.09 


gold 


SN 2003az 


1.265 


45.20 


0.20 


0.25 


gold 


SN 2003ak 


1.551 


45.30 


0.22 


0.86 


gold 


SN 2003bd 


0.67 


43.19 


0.28 


0.27 


gold 


SN 2003be 


0.64 


43.07 


0.21 


0.23 


gold 


SN 2003dy 


1.340 


45.05 


0.25 


0.54 


gold 


SN 2003es 


0.954 


44.28 


0.31 


0.07 


gold 


SN 2003eq 


0.839 


43.86 


0.22 


0.22 


gold 


SN 2003eb 


0.899 


43.64 


0.25 


0.26 


gold 


SN 20031v 


0.94 


43.87 


0.20 


0.15 


gold 



*Peculiar velocity of 400 km s ^ included for all SNe and 2500 km s ^ if ^ from SN. 
"Distance normalization is arbitrary; see Appendix 
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sample are shown as diamonds, //.ST-discovered SNe la are shown as filled symbols. Bottom 
panel: weighted averages in fixed redshift bins are given for illustrative purposes only. Data 
and models are shown relative to an empty Universe model (0 = 0). The fit statistics 
for each model are listed in Table 4. 

Figure 8: Joint confidence intervals for (Qm^^a) from SNe la. The solid contours 
are results from the gold sample of 157 SNe la presented here. The dotted contours are 
the results from Riess ct al. (1998) illustrating the earlier evidence for f^A > 0. Regions 
representing specific cosmological scenarios are illustrated. Contours are closed by their 
intersection with the line Qm = 0. 

Figure 9: Joint confidence intervals for Qm and a static equation of state for dark 
energy, w. In the left-hand panel, constraints from the gold SN la sample (dotted contours) 
are combined with a prior of Qm = 0.27 ±0.04 to yield the solid contours. In the right-hand 
panel, the same SN constraints are combined with those from WMAPext and 2dfGRS to 
yield the solid contours. 

Figure 10: Joint confidence intervals derived from SN samples for a two-parameter 
model of the equation of state of dark energy, w{z) = Wo + w'z. For each panel, constraints 
from a SN sample is combined with the prior, Qm = 0.27 ± 0.04, to yield the indicated 
confidence intervals. The position of a cosmological constant (—1, 0) is indicated as a filled 
symbol. The lower-right panel shows the impact of adding or subtracting a systematic error 
in distance modulus of 0.05^ mag to the gold sample. 

Figure 11: Comparison of composite rest- frame light curves of SNe la. Each SN la is 
individually transformed to the rest frame (K-corrected; corrected ior 1 + z time dilation). 
Each is normalized by subtraction of the peak magnitude and its date. Data from SNe la 
with 2; < 0.1 are shown as open symbols; data from the iJST-discovered SNe la are shown 
as filled symbols. 

Figure 12: Light-curve shape distributions of low- redshift and high-redshift SNe la. 
Histograms of the MLCS2k2 parameter A (relative peak visual luminosity) are shown for 
SNe la with 2; < 0.1 and the if^T-discovered SNe la. 

Figure 13: Rest-frame Bmax — Vmax and Umax — B^ax colors of SNe la versus redshift. 
Expected, unreddened colors of SNe la are shown as dotted lines. 

Figure 14: Histograms of the color data shown in Figure 13. The distributions from 
three different redshift bins are shown as indicated. 

Figure 15: Comparison of individual distance difference estimated by the MLCS2k2 
and BATM methods for the HST and ground-discovered SNe la. The zero-points of both 
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methods are normalized by using the same set of SNe la. 

Figure 16: Predicted lensing magnifications of SNe la and of random positions in 
the CDF-S and HDF-N. For the SNe la discovered in the GOODS fields, the expected 
magnification was calculated using a multiple-lens plane formalism, with estimates 
of foreground lens redshifts and masses derived from the GOODS catalog. Expected 
magnifications were also calculated for 100 randomly selected positions (redshift and 
angular position). The solid and dotted lines show redshift bin averages and dispersion, 
respectively. SNe la found in the GOODS survey and in other HST searches are indicated. 

Figure 17: Correlation of the predicted magnification and the best-fit cosmological 
model residual for individual SNe la. The predicted magnifications are as described in 
Figure 16. The residuals arc the difference in distance modulus as predicted from the 
best-fit model {Q.m = 0.3, Qa = 0.7) and as observed. The empirical correlation is expected 
to be unity (if the lens light traces their mass) and is shown for the whole sample and 
without SN 1997ff, the SN with the largest residual and predicted magnification. 



-51 - 



References 

Aguirrc, A. N. 1999a, ApJ, 512, L19 

. 1999b, ApJ, 525, 583 

Alard, C, & Lupton, R. H. 1998, ApJ, 503, 325 
Barris, B., et al. 2004, ApJ, in press 

Bemtez, N., Riess, A., Nugent, P., Dickinson, M., Chornock, R., Filippenko, A. V. 2002, 

ApJ, 577, LI 
Bennett, C, et al, 2003, ApJS, 148, 1 

Blakeslee, J. P., et al. 2003, ApJ, in press (astro-ph/0302402) 
Boehm, A. et al. 2003, A&A in press, (astro-pli/0309263) 

Budavari, T., Szalay, A. S., Connolly, A. J., Csabai, I., & Dickinson, M. 2000, AJ, 120, 1588 

Caldwell, R. R., Dave, R., & Stcinhardt, P. J. 1998, Ap&SS, 261, 303 

Caldwell, R. R., Kamionkowski, M., & Weinberg, N. N. 2003, Phys. Rev. Let., 91, 71301 

Carroll, S. M., Press, W. H., & Turner, E. L. 1992, ARA&A, 30, 499 

Chapman, S. C, Blain, A. W., Ivison, R. J., & Small, 1. 2003, Nature, 422, 695 

Clochiatti, A., et al. 2004, in preparation 

Coil, A. L., et al. 2000, ApJ, 544, Llll 

Colgate, S. 1979, ApJ, 232, 404 

Cowie, L.L., Barger, A. J., Hu, E. M., Capak, P., & Songaila, A., 2004, AJ, submitted 

(astro-ph/0401354) 
Deffayet, C, Dvali, C, & Gabadadze, G. 2002, Phys. Rev. D, 65044023 
Di Pietro, E., & Claeskens, J. 2003, MNRAS, 341, 1299 
Drell, P. S., Loredo, T. J., & Wasserman, I. 2000, ApJ, 530, 593 
Falco, E., et al. 1999, ApJ, 523, 617 

Farrah, D., Fox, M., Rowan- Robinson, M., & Clements, D., 2004, ApJ, in press 
Fernandez-Soto, A., Lanzetta, K. M., & Yahil, A. 1999, ApJ, 513, 34 
Filippenko, A. V. 1997, ARA&A, 35, 309 
. 2001, PASP, 113, 1441 

. 2003, in From Twilight to Highlight: The Physics of Supernovae, ed. W. Hillebrandt 

& B. Leibundgut (Berlin: Springer- Ver lag) , 171. 
. 2004, in Carnegie Observatories Astrophysics Series, Vol. 2: Measuring and Modeling 

the Universe, ed. W. L. Freedman (Cambridge: Cambridge Univ. Press), in press 

(astro-ph/0307139). 

Filippenko, A. V., Li, W. , Treffers, R. R., & Modjaz, M. 2001, in Small-Telescope 
Astronomy on Global Scales, ed. W. P. Chen, C. Lemme, & B. Paczynski (San 
Francisco: ASP), 121 

Filippenko, A. V., & Riess, A. G. 2001, in Particle Physics and Cosmology: Second Tropical 
Workshop, ed. J. F. Nieves (New York: AIP), 227 



-52- 



Preedman, W., & Turner, M. 2003, Rev. Mod. Phys. CoUoquia, in press (astro-ph 

/0308418) 

Fnichter, A., & Hook, R. 1997, SPIE, 3164, 120P 
Garnavich, P. M., et al. 1998, ApJ, 509, 74 
. 2002, BAAS, 78.09 

Giavalisco, M., et al. 2003, ApJL, in press (astro-ph/0309105) 

Gilliland, R. L., Nugent, P. E., & Phillips, M. M. 1999, ApJ, 521, 30 

Gilliland, R. L., & Riess, A. G. 2002, in HST Calibration Workshop, ed. S. Arribas, A. 

Koekemoer, & B. 
Goobar, A., Bergstrom, L., & Mortsell, E. 2002, A&A, 384, 1 
Gunnarsson, C. 2004, JCAP, submitted (astro-ph/0311380) 
Gwynn, S. D. J., & Hartwick, F. D. A. 1996, ApJ, 468, L77 

. 1996a, AJ, 112, 2391 

. 1996b, AJ, 112, 2398 

Hamuy, M., Trager, S. C, Pinto, P. A., Phillips, M. M., Schommer, R. A., Ivanov, V., & 

Suntzeff, N. B. 2000, AJ, 120, 1479 
Hatano, K., Branch, D., & Deaton, J. 1998, ApJ, 502, 177 
Hoflich, P., Wheeler, J. C., & Thielemann, F. K. 1998, ApJ, 495, 617 
Holz, D. E. 1998, ApJ, 506, LI 
Jha, S. 2002, Ph.D. thesis. Harvard University 
Jha, S., Riess, A. G., & Kirshncr, R. P. 2004a, in preparation 
Jha, S., et al. 2004b, in preparation 
Jha, S., et al. 2004c, in preparation 
Kallosh, R., & Linde, A. 2003, JCAP, 302, 2 

Kallosh, R., Kratochvil, J., Linde, A., Linder, E. V., & Shmakova, M. 2003, JCAP, 310, 15 

Kantowski, R. 1998, ApJ, 507, 483 

Knop, R., et al. 2003, ApJ, in press (astro-ph/0309368) 

Kochanek, C. S., 1996, ApJ, 473, 595 

Koekemoer, A., et al. 2004, in preparation 

Leibundgut, B. et al. 2004, in preparation 

Leibundgut, B. 2001, ARAA, 39, 67 

Lewis, G. F., & Ibata, R. A. 2001, MNRAS, submitted (astro-ph/0104254) 

Linder, E. V. 2003, Phys. Rev. Lett., 90, 91301 

Linder, E. V., & Huterer, D. 2003, Phys. Rev., D67, 81303 

McLean, L S., et al. 1998, Proc. SPIE, 3354, 566 

Metcalf, R. B., & Silk, J. 1999, ApJ, 519, 1 

MortseU, E., Gunnarsson, C., & Goobar, A, 2001, ApJ, 561, 106 (MGG2001) 
Ninomiya, Yaqoob, T. & Khan 2003, private communication 



-53- 



Nonino, M. et al., 2004, in preparation 

Nugent, P., Kim, A., & Pcrlmuttcr, S. 2002, PASP, 114, 803 

Okc, J. B., & Sandage, A. 1968, ApJ, 154, 21 

Oke, J. B., et al. 1995, PASP, 107, 375 

Paerels, F., Petric, A., Telis, G., & Helfand, D. J. 2002, BAAS, 201, 9703 
Pain, R., et al., 2002, BAAS, 43.02 

Parker, L., Komp, W., & Vanzella, D. 2003, ApJ, in press 
Parker, L., & Raval, A. 1999, Phys. Rev., D60, 123502 
Peebles, P. J., & Ratra, B. 2003, Rev. Mod. Phys., 75, 559 

Perlmutter, S., Sz Schmidt, B. P. 2003, in Supernovae & Gamma-Ray Bursts, ed. K. Weiler 

(New York: Springer, Lecture Notes in Physics), in press 
Perlmutter, S., et al. 1997, ApJ, 483, 565 
. 1999, ApJ, 517, 565 

Poznanski, D., Gal- Yam, A., Maoz, D., Filippenko, A. V., Leonard, D. C., & Matheson, T. 

2002, PASP, 114, 833 
Rana, N. C. 1979, Ap&SS, 66, 173 
Rana, N. C. 1980, Ap&SS, 71, 123 
Riess, A. G. 2000, PASP, 112, 1284 
Riess, A. G., et al. 1998b, AJ, 116, 1009 

. 1999a, AJ, 117, 707 

. 1999b, AJ, 118, 2675 

. 2001, ApJ, 560, 49 

. 2003, ApJ, 600, in press 

Riess, A. G., Nugent, P. E., Fihppenko, A. V., Kirshner, R. P., & Perlmutter, S. 1998a, 
ApJ, 504, 935 

Riess, A. G., Press, W. H., & Kirshner, R. P. 1995, ApJ, 438, L17 

. 1996a, ApJ, 473, 88 

. 1996b ApJ, 473, 588. 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Schmidt, B. P., et al. 1998, ApJ, 507, 46 

Scranton, R., et al. 2003, Phys. Rev. Lett., submitted (astro-ph/0307335) 

Sirianni, M. et al. 2004, in preparation 

Smith, R. C., et al., 2002, BAAS, 78.08 

Spergel, D. N., et al. 2003, ApJS, 148, 175 

Strolger, L.-G., et al. 2004, in preparation 

Sullivan, M., et al. 2003, MNRAS, 340, 1057 

Suntzeff, N., et al. 2004, in preparation. 

Tammann, G. A. 1979, in ESA/ESO Workshop on Astronomical Uses of the Space 



-54- 



Telescope, ed. D. Macchetto, F. Pacini, & M. Tarenghi (Geneva: ESQ), 329 
Tonry, J. T., et al. 2003, ApJ, 594, 1 

Treu, T., Stavelli, M., Cascrtano, S., Moller, P., & Bertin, G., 2002, ApJ, 564, 13 
Turner, M., & Riess, A. G. 2001, ApJ, 569, 18 
Visser, M. 2003, gr-qc/0309109 

Wambsganss, J., Gen, R., & Ostriker, J. P. 1998, ApJ, 494, 29 
Wirth, G. D., et al., 2004, AJ, submitted (astro-ph/0401353) 
Wright, E. L. 2002, BAAS, 161.17 (astro-ph/0201196) 
Ziegler, et al. 2002, ApJ, 564, 69 



-55- 



Fig. 1. — Figure 1. See attached jpeg image of SN hosts. 
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